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Images* A few subtle but mappahle thermal-inertia anomalies coincide with areas of 
helium leakage believed to be associated with deep oil and gas concentrations. 
Similar changes were not found In areas of uranium deposits possibly due to theAt 
minor and dlacontlnuout nature at satellite resolution. 

The most Important results Involved delineation of tectonic framework elements 
some of which were not previously recognlstd. Thermal and thermal-inertia Images 
also permit mapping of gcomorphlc testursl domains. A thermal lineament not on 
siiatlng geologic maps or detected on Landsst Images appears to reveal a basement 
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hydrothermal faatures could not be distinguished with any confidence. 
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Inertia mapping algorithm, a fast and accurate Image-registration technique, and an 
afficlent topographic-slope and elevation-correction method. 
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1.0 PREFACE 


1.1 OBJECTIVES 

The principal objective of thii study was to investigate applications of 
HCMM (Heat Capacity Happing Mission) satellite data in detecting and caapping 
geologic features for energy~resource and mineral-'deposit studies. Other, 
related objectives involved the development of new techniques and approaches 
in thermal modeling and image processing. 

1.2 SCOPE 

The analysis was somewhat restricted by the limited number of sequential 
day/night image pairs free of major atmospheric/weather problems. For the 
Powder River Basin, Wyo. area, a single thermal” inertia image was formed using 
20 August 1978 data. Four additional nighttime scenes were used to examine 
geologic formation boundac.es and thermal lineaments. The analysis of the 
Cabeza Prleta, Ariz., area was done using a thermal-inertia image constructed 
from data acquired April 3 and 4, 1979. No successful U-2 alrcraf data 
acquisition flights were conducted over these two sites so that comparison of 
different resolution thermal data was only conducted using USGS aircraft 
data. Despite all thes« limitations, significant geologic information was 
derived in this study, and the results suggest the Importance of a follow-on 
thermal satellite experiment for Improved mineral and energy resource explora- 
tion. 
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1.3 CONCLUSIONS 


Octplte Halted d«ta» InvastlgatlonB In the Powder River Basin area of 
•astern Wyonlng end adjacent States clearl showed that geologic units sk 
narrow as two or three resolution eleraents, but of Boderate to high thermal- 
inertia contrast against surroundings, can be discriminated in optimal 
Iraagest It appears likely that subtle facies differences in sedimentary 
basin-’flu units can be delineated and mapped using satellite thermal-inertia 
Images, especially If sequential Images can be obtained during a drying cycle 
after rain or snoj. A few subtle but mappable thermal-inertia anomalies 
coincide with areas of anomalous helium In soil gas believed to Indicate 
leakage from deep oil and gas concentrations; the presence of thermal-inertia 
anomalies suggests that gas leakage has produced chemical changes and cementa- 
tion at the surface. Such changes also are known to be associated with shal- 
low uranium deposits and changes were looked for but not found in the thermal- 
inertia images; It Is thought that the surface changes In this area are too 
minor and discontinuous to be detected from satellite. 

The most consistently practical and important results involve delineation 
of tectonic framework elements such as lineaments bounding apparent structural 
blocks. These commonly can be seen even in less-than-optimal data. One pair 
of major thermal lineaments In the southern Powder River Basin seems to define 
structures not previously recognized but consistent with, and adding impor- 
tantly to, an emerging atory of basement-block movements and their direct 
Influence on sedimentation, which in part controls the occurrence of large oil 
and gas resources. One of these lineaments matches up with aeromagnetic nup 
data and appears to reveal a basement dlscontinuivy which underlies the famous 
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HoBtitakc Mine In the Black Hille and a lone of Tertiary Igneous activity. 
Along with the newly Identified llneanents» the thermal Images also permit 
aapping of geomorphlc textural domains. The ge^ ' ^>£10 significance of these Is 
not yet understood, but it seeni likely tt they connote structural and 
lithologic conditions which affect or control local ground-water regimes. 

Similar applications of HCMM data to the Cabeza Prieta, Arlz., area 
illustrate the potential of using themal-lnertia data for discrimlnac ion 
between extrusive and Intrusive rocks and for detecting differences in the 
mafic content of volcanlcs. Other results included detection of differences 
among surflclal units - tentatively ascribed to changes in soil-moisture 
retention, dls:overy of discrepancies in existing geologic maps, and possible 
application of the thermal-inertia technique to mapping burled pediments. 

Extension from beyond the originally proposed study areas to Yellowstone 
National Park was made to examine the usefulness of HCMM data in geothermal 
studies. Although we found that the night-thermal data could not be used with 
any confidence to distinguish surface hydrothermal features, we did detect 
additional structural Information concerning the outline of the caldera which 
Is the source of the volcanic heat. This reinforces our conclusion that a 
raajov utility of these data Is in providing information about local-area or 
regional tectonic framework. 

We have also made significant advances In modeling analysis and image- 
registration techniques. A thermal-inertia mapping algorithm has been devel- 
oped based on a new method to derive the regional meteorologlc parameters 
solely from tht satellite data. An algorithm for dete.-rolnlng the sensible- 
heat flux from ground-station data was also constructed. Simple forms for 
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four of the atBOsphorlc flux ttrns war* eonstructod froa flold waiurtnentR 
■ode during circuattancas whan aatallita data ara llkaly to ba aoat useful. 
Thasa forms aliainate the need for extensive continuous ground station data. 
Also, a aathod to correct thermal and thermal-inertia data for elevation 
variations In sky and aolar flux was deterstlned. In addition^ we have devised 
a fast topographic adjustment algorithm which can be used Ivk conjunction with 
digital terrain data to correct the thermal-inert la Image for simple topo- 
graphic slope affects. Flnallyi a fast Image registration technique was 
developed that proved to be considerably more accurate than the NASA regis- 
tered products. 

Our analysis of the HCMM data has resulted In the recognition of features 
which suggest tha existence of previously unmapped and unknown geologic struc- 
tures. Their relationship to other geophysical and geochemical data provides 
Important information for a basic resource-exploration strategy. Addition- 
ally, substantial progress has been made in modeling and Image-processing 
techniques. This report covers new areas and represents significant advances 
In the processing and Interpretation of thermal satellite data and in the 
integration of thermal-infrared data in regional ge<ilogic exploration. 

1.4 RECOMMENDATIONS 

From our experience to date, we would recommend that serious consider- 
ation be given to a follow-on thermal satellite mission with these general 
characteristics. 

1. The current NEAT of HCMM seems adequate for most regional studies. 
Higher thermal resolution does not appear necessary. 
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2. Sobs Incrssto in ground rotolutlon (possibly 100>200 n) would bo 
useful: however, there ere trade-offs to consider here* The 500-n resolution 
froa HCMN has proven very useful for regional structures - it does not appear 
proaising for detecting alteration* 

3. Some increase in the repeat times over a site is desirable* The HCMM 
data we have seen have often baffled us because of changing oeteorologic 
effects* The increased repeat time would enhance the chances of "stable'* 
clear" conditions and also provide coverage of regions under several meteoro- 
logical and soil moisture conditions* The repeat time involves the orbit 
parameter selection; a 5-10 day repeat of coverage would be desirable* 

4* The current overflight times of HCMM appear appropriate for geologic 
analysis* It should also be noted that the daytime maximum represents an 
optimum time to acquire multispectral thermal measurements as well* 

5* Our analysis of HCMM data requires registration of day and night 
images and subsequent registration to a topographic base* The registered data 
provided by NASA often contained large registration errors* An essential 
requirement for analysis of these data is that the clear scene images be 
registered (day/night images) to a pixel, and to digital terrain data* If 
this registration accuracy cannot be achieved routinely, it is recommended 
that registered products not be provided to users* 
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2.0 INTRODUCTION 

Th« UCKM dat« vhleh wt havt MMlntd have providtd ut vlth uniqut gtolog'* 
it InforMClon which li both eospltx to analyct and difficult to explain. In 
•oac caati wt art astoniahad to aaa aubtla diatinctiona of atructura and 
gaologic Mtariala that ara not found on datailad gaologlc laapa. In other 
caaaa wa ara unabla to dif farantiata wldaly diaaiailar gaologlc Mtariala or 
idantlfy faaturaa which ara claarly daacribad on ragional gaologlc naps or 
Landaat iaagaa. Coaiaonlyi geologic faaturaa ara claarly diaplayad on a aingle 
loaga or part of an iaaga and not on othara. Further complicating the analy- 
aia of thaae data haa been the axperimantal nature of the aatallita miaaion, 
which haa introd'^cad both aignificant time laga batween data acqulaltion and 
interpretaMon, and unique conatrainta in image raglatration and data calibra- 
tion. 

Thia report compriaea pieces to a purzle— *a puzzle with tantalizing new 
information tut sufficient gape to preclude a complete overall asseasment. We 
have examined the geologic implications of the HCMM data in the Powder River 
Basin of Wyoming and subsequently extended the analysis to the Cabeza Prieta 
area in Arizona. Enhanced nighttime images, thermal-inertia images based on 
our new algorithm, and several profiles of various data across the basin are 
presented in our analysis and subsequently compared with other map data (geo- 
logic, Mgnetic field, known areas of oil, gas, and uranium occurrence, helium 
anomalies, and ground water). Of equal importance with the geologic interpre- 
tation haa been progress in modeling. We provide several new algorithms: 
thermal- inertia Mpping, estimation of regional meteorological Information 
from the fundamental remote sensing data, registration of satellite day-night 
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tltvtclon correction of thtnul and thtnMl-inortla data, and dctartsi- 
nation of aontlbla-haat flux* The concludlnt aactlon of thla report addret«ef 
our iMge proeeerl'^ technique! together vith a Hating of the conputer pro- 
Craai used. 

2.1 GEOLOGIC SETTING 

The atudy areas arc the Powder River Basin and environs In eastern Wyo- 
■ing'Montsna and the adjacent Dakotas and the Cabesa Prlcts Range in southwest 
Arisons. The Powder River Basin (1st. 42°-45° N. , long. 103°-107° W.) has 
large potential for coal, oil and gas, and uranium, and accordingly is now the 
target of several major geologic, water-resourcc, and land-use mapping pro- 
jects. Covering an area of about 250 x 400 km. It is a semi-arid region of 
rolling low hills typically with thin to moderate grass and sage cover. 
Tertiary rock units (Port Union and Wasatch Formations) in the central arcs of 
the basin, vherc the energy resources are known to occur, are exposed on 
scales sufficient for ''stellite measurements. The lower part of the Fort 
Union is sandstone exposed in belts 4 to 10 km wide; the upper Fort Union part 
is siltstonc with major coal beds exposed in belts 10 to 30 km wide; and the 
Wasatch is siltstone and claystone covering areas 30 to 60 km wide. 

The Cabeza Prleta Range in Arizona (lat. 32°-33^ N. , long. 112°-115° W.) 
is s proposed Wilderness area, and the USGS has begun a program to define the 
g<eology and miner-l-resource potential of this virtually unmapped area. The 
State geologic map shows the area to contain granite, schist, mafic volcanic 
rocks, and alluvium. It lies very near the major mineral district at Ajo, 
Ariz., and contains old prospect developments in hydrothermally altered 
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ground. itcauio tht «r«a has bttn withdrawn from public accais aince World 
War 11, tha gaology and aineral potantial art bartly known and thua the area 
it contidtrtd rtlativaly iapritant for aodern study. 


2.2. APPROACH 

Initial data interpretation was performed by visual pattern recognition 
of areas significantly different than their surroundings. Several 
1 : 1,000,000-scale photographic enlargements were made of the NASA thermal and 
reflectance images to match them to many of the other existing data map prod- 
ucts (geologic, geophysical, topographic, Landsat lineaments, and so on). 
Because of the X-Y distortion and the large magnification, this method was not 
entirely satisfactory for detailed study. We then generated film products 
from the computer-compatible tapes (CCT's) and used a zoom transferscope with 
X-Y stretch capability to register the projected images onto a stream-network, 
map. In mont local areas, this permitted plotting of features to within one 
to three pixels of their true ground position. Few of the interpretations 
based on the temperature boundaries or on other features are significantly 
affected by this degree of mlslocation. 

During this stage of the investigation, we also discovered that the NASA- 
supplied AT and thermal-inertia images contained artifacts, such as double 
drainage, indicating misregistration in parts of these images of several to 
many pixels. Consequently we developed a registration algorithm (Watson and 
others, 1981b) which registers data to within two pixels. Also as part of the 
Investigatlrri, we developed a new thermal- inertia algorithm (Watson, 1981a) 
which was employed during the remainder of our analysis. Correlation of these 
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ragittarad laaga prodMCta with tha othar gaologlc-gaophyaical data was per- 
foraad prlaarily by using optically anlargad projactions of tha data on digi- 
tally anlargad iaaga product at a scale of approxluataly li200,00U. This 
anablad us to oxaoine aubtla faaturts at tha plxtl level and also to enploy 
our full profiling and histogram capability at the full dynamic range and 
resolution of the digital isuige data. We have thus been able to quantify many 
of the scene differences which had been observed on various image products. 

We also examined the use of color-coded images for enhancing subtleties 
in the scene contrast. Cenerally this provided a more obvious demonstration 
of differences but did not appear to add any new information. Toward the end 
of our study, howe.rer, a color-coded thensal-inertia image was produced which 
provided a new - if unexplained “ perspective of the scene. A north- 
noiineist’^trending rectangular pattern of ground, surrounding the Black Hills 
and roughly corresponding with major changes in the drainages of the Yellow- 
Btone and Hlssouri Rivers was observed In retrospect, this feature can now 
also be seen on the black and white products. The color-coding of the image 
was also a very useful tool for quickly determining the numerical range of 
values. With an appropriate color-scale and u^Ung a high-powered magnifica- 
tion lens, it was possible to determine the thermal-inertia ranges of many 
geologic units quickly. In both these respects the color-coding can be re- 
garded as a useful but not csrential element in the analysis. 
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2. 3 RESULTS 


2* 3*1 ThT«al"intrtl« — pptrg for dltcrlwlnatlon of gtologic futures 

2* ^.1.1 Dtlincatlon and tubdlvlilon of geologic units 

Discrimination studies have been somewhat limited by the lack of sequen* 
tlal day~nlght image pairs free of major atmospheric/weather problems, but the 
20 August 1978 set of day (AAOl 16-20010-1, 2 and AA0116-20020-l,2) and night 
(AAO 11 6-09040-3) scenes and our constructed temperature-difference (<^1} and 
thermal- inert la images show several geologically significant features. Other 
good nighttime data from 30 July 1978 (AA0095-09 170-3), 5 September 1978 
(AAOl 32-09050-3), 27 September 1978 (AAOl 54-091 90-3), and 10 June 1979 
(AA0410-08450-3) passes have been used for delineation of several geologic 
bov’Tidarles, including some not within the area of, or not identified in, the 
20 August 1978 data set* 

A measure of discrimination capability, using optimal images ana selec- 
tion of geologic units which contrast well with their surroundings, is the 
clear delineation of the Mesaverde Formation* South of the Bighorn Mountains, 
In the vicinity of the towns of Midwest and Edgerton Wyo* , the Mesaverde - a 
relatively massive sandstone — crops out between the Fox Hills Sandstone and 
the Cody Shale* On the night image (fig* 1) for 30 July 1978, all three units 
are relatively warm; however, the Mesaverde can be traced as a distinctly 
warmer unit (l/4°-l/2° C), along at least 40 km of strike length in which a 
fold nose is clearly defined (fig* 2)* The Mesaverde outcrop here is 1-2 km 
wide or 2-4 resolution elements (pixels). 
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Southwest of that area, south of ths town of Powder River, within a large 
araa of Cody Shale, is a aandstone inlier of Frontier Foraation surrounding 
older shale and sandstone of the Clovcrly and Morrison Forvations* In the 3U 
July iaage, the Frontier is clearly a warm annulus around a cool center of the 
other units (fig* 1 and 3)* The units in the cool central area are about b km 
wide and the werm Frontier annulus is 3.5 ka wide* What is additionally 
interesting in this area is the apparently clear definition of a similar but 
smaller such feature to the northwest which does not match the shape of the 
contacts on the siost recent geologic map (Love and others, compilers, 1935). 
This feature has not been field checked* 

Another measure of discrimination is found in the area of the Pumpkin 
Buttes* These are very sharply defined in the 3 September 1978 night image 
(fig* 4). North Butte is about 3 km vide, and the combined topographlc/geo- 
logic prominence of Middle and South Buttes measures about 4 by 8 km* These 
are warmer than their surroundings, aa is expected of tuffaceous sandstones of 
the Whiti River Formation, dense and resistant enough to form buttes where 
erosional remnants lie upon the softer sandstones and mudstones of the Wasatch 
Formation* Once again, when the geologic map was projected onto this HCMM 
scene, differences were observed. Unfortunately, on our only thermal- inert la 
image, clouds were present over the Buttes preventing the observation of the 
expected thermal-inertia contrast. 

Definition of geologic features is highly variable from pass to pass and 
within single passes. The 30 July image (fig. 1) is excellent for the north 
half of the Powder River Basin and the areas west and southwest of the 
Basin. The image appeara virtually washed out in the south half of the Basin, 
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showing high tomptrscurss In sn srss conspicuously cool in all the other 
iaagest Hscausc t.ha image was examined long after it was acquired, it was not 
possible to obtain detailed field meteorological data. This is a generic 
problem in dealing with transient and local phenomena that commonly affect 
thermal surveys. Air temperature and precipitation data rhow fairly similar 
conditions at the 28 weather stations throughout the scene. From these data, 
the intrascene variations cannot readily be ascribed to local weather/moisture 
changes; however, the NOAA and hMSP (Defense Meteorological Satellite Program) 
Satellite data show that a major weather front had recently passed through the 
basin. Such intrascene differences are less pronounced or even absent in 
other passes, but none of the others expresses quite the same degree of geo- 
logic feature definition as the good portions of the 30 July image. 

The multiple data sets of 20 August contain considerable geologic infor- 
mation, especially in comparing patterns seen variously in the day thermal, 
night thermal, AT, and thermal-inertia images. The day thermal image (fig. 3) 
shows large areas of warm ground north and east of the Black Hills. These do 
not correspond to lithologic subdivisions on any available geologic maps, nor 
to any patterns of weather across the scene during the previous few days. 
Small individual features of interest in the image are cool areas around the 
Tongue River and in a belt of small patches trending north-south up the center 
of the southern part of the basin. The very warm drainage area west of the 
Black Hills and the warm area south of the Black Hills are also noteworthy. 
The night image (fig. 6) offers busier patterns of finer scale definition, 
dominantly related to the topographic character of local areas- Much, but not 
all of the high ground between streams, is conspicuously warm. The long, 
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■tr«l|ht Fowdcr Rlv«r Mparatti Man gv^ound on let o«at flank fro* cooltr 
ground on ita waat flank. Thia la not juat an affact of aaat- va. wtai -facing 
alopta» bacauaa aoat othar atraan araaa do not ahow tha aant affact. The 
north-aouth bait In tha aoutharn part of tha iaain ahowa aa vary warn patchea 
of ground. Vary warn araaa ring tha aaat flank of tha Black Hllla and occur 
Juat to tha north throughout tha Baar Lodga Mountaina, but thaaa do not cor re*- 
apond to tha nappad gaology. Othar eontraat atratchaa vara triad including 
color alicing but tha corralationa of tanparatura araaa with gaology did not 
improve. What doaa appaar to be trua» however, ia that within the Powder 
River Baain, conaplcuoualy warm araaa arc much more abundant in the north 
half. For the moat part, thia ia a reault of gn^atar diaaactlon of the ter- 
rain and more axpoaura of bedrock, aa compared with few outcropa and abundant 
windblown eand veneer in the aouth half. 

Analyaic of the thermal-inertia imagaa, derived from our regiatration and 
modeling algorithms, ahowed that the Tongue River areas of cold ground in 
daytime are, in fact, areac of high thermal inertia (2000 Thermal inertia 
units (TIU); 1 TIU ■ I W aec^^^ m”^). These areas (fig. 7) correspond quite 
well to areas mapped by Raine* (Raines and others, 1978), using computer 
enhancements of Landsat images. They were mapped as the coarsest, sandiest 
lithofacies unit in the baain, which are relatively indurated and resistant 
and should crop out beat and, depending on moisture conditions, should have 
the highest thermal inertia of the subunits in the Wasatch and Fort Union 
Formations. Other auch correlations exist for several areas of this facies 
aoutheastward toward the Black Hills. Areas of Che Wasatch and Fort Union, 
sampled from various parts of the Basin and which appear to be representative 
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of larttr ■urroundini artai h«vt th«nul*in«rtlas of 1525 TIU and 1A50 TIU, 
raapactivoly Howovtr» tha thonMl*lnartia of those two units varies consid** 
•rably throughout tho basin and in somo cases values for these two units are 
•tatistically inseparable. A possible explanation is that these units genvr- 
ally forn low or flat topography where windblown sand obscures the underlying 
geology in extrenely irregular (and unmapped) patterns. Another reason is 
that these units retain moisture differently and longer than sandier facies, 
and thus may have great irregularities in both thermsl'inertia values and wind 
cooling patterns. If they are slightly wet, and not cooled by surface winds, 
their thermal inertia will be higher. 

The north-south belt of high thermal inertia (fig. 7) was initially 
thought to correspond to burned ground over ancient natural coal fires. ITie 
night thermal image (fig. 4) was carefully registered to a base map and a 
composite map was made showing the areas of clinkers (as determined from a 
color ratio composite Landsat image) and the warm areas on the HCMM image 
(fig. 8a). We then examined the thermal-inertia image and determined that the 
clinker areas in fact have an intermediate thermal inertia (1300 Tl'J) and the 
N-S belt of warm ground in the night image Just east of the clinker hills has 
a higher thermal inertia (1500 TIU). This north-south belt has been mapped in 
detail and the surface ology f/rovides no clue as to why these areas have 
high thermal inertia. This docs not conform with conditions produced where 
windblown sand accumulates in the lee of topographic highs. It is suspected 
that the highly fractured clinker hills are readily drained of their near- 
surface moisture and this ground water tends to pond Just eastward in the 
direction of tvormal drainage, causing an increase in thermal inertia. This 
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hypothaiii 1m glvtn cradancc by an axamlnation of the NURE gaama-ray profiles 
In thlN area. Thaae araaa coincide with lowe in the total gavaa-ray measuti"- 
■anti (fig. 6b) aa would be axpacted for araai of higher noiatura content. 

A rough order of aagnitude aatiouite from theae data la that the anomalous 
araaa are aaaoclated with a 20 percent increaxe in thermal inertia amt a S ii< 
10 percent decraaae in the total count valuea. The tltermal inertia of koIIk 
increaaea rapidly with incrwaalng aoil nolature content and the effect can he 
aatimated for low moiature content! by conaidering only that increase due to 
denaity and apeclfic heat capacity. The ratio of the fractional change in 
thermal inertia to density la Just one half the ratio of the apecific heat 
capacity of water to soil or approximately 2.3. Thus a 20 percent increasi* in 
thermal inertia could be produced by a aoil moisture change which increahes 
the density by 8 percent (and decreases tite total count by an equivalent 
percent ). 

To examine the basin further, two northwest-southeast profiles across tl>e 
20 August image (fig. b, profiles A-A\ B-U') were constructed. These pro- 
files enabled us to look in detail (pixel level) at variations in tlieniinl- 
Inertia values and to examine relationships between temperature or thermal- 
inertta patterns and topography. Topographic data were taken from 1:230,00;' 
USOS base maps with a contour interval of 200 feet. Figures 9 through I show 
profiles of thermal-inertia , elevation, and topographic gradient along linos 
A-A' and B-B*. 

The profiles on A-A' have several Interesting features. The line begins 
In the Wasatch Formation at the northwest end, and thermal-inertia vaUns 
(fig. 9) decline Into a broad low, about coincident with the Powder River 
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Ffgurs 8b.-*** Comptrfaon batwaan profllaa of total 

oount gamma ray and tha anomaloualy warm 
araaa.Tha atralght Hnaa Indloata tha 
gaographio poaltlon of tha profllaa. 







drainage aeen ao clearly in the topographic profile (fig* 10, 11), and then 
riae. This low approxinuitely aarks the basin axis, and the adjacent slopes of 
the thet'mal”lnertia profile show the character of the upper part of the 
Wasatch on either side of that axis. The topography itself is somewhat dif- 
ferent on opposite flanks of the Powder River drainage, and this probably 
explains the previous observation that the night- temperature image showed the 
two flanks differently, even though the thermal-inertia image indicated the 
two flanks to be underlain by similar material* The flanks have different 
slopes and bed dips, and the western flank generally is dissected more sharply 
and deeply than the eastern flank. At the next large drainage east of the 
Powder River, the thermal-inertia profile breaks sharply, suggesting either a 
previously unmapped lower unit of the Wasatch or a sharp change to the some- 
what finer grained facies which has been noted in the lower part of the forma- 
tion. Along this profile the Fort Union Formation has a roughly estimated 
average value of 1423 TIU, as compared with an equally rough, general average 
of the Wasatch of 1623 TIU. This difference is about what would be expected 
from the compositions of the two formations, although they are rather nonuni- 
forra on the scale of the whole basin. The "typical" areas of Fort Union and 
Wasatch that were sampled gave values of 1430 TIU and 1323 TIU, respective- 
ly. A sharp break in the thermal-inertia profile occurs between the two 
formations, but it falls 4 to 3 km west of the contact as shown on the geo- 
logic map. A break or dip also occurs in the profile at the contact of the 
upper (Lebo Shale) and lower (Tullock) members of the Fort Union, but overall 
the members have about the same thermal Inertia. The elevation profile shows 
a marked change in character of topography from Wasatch to Fort Union, as does 
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th« topographic-gradlant profile* The wera north-eouth tone le e nerrow but 
proalncnt therael'lnertle epike in the lower pert of the Lebo Shele Member of 
the Fort Union. The Pierre Shele hee the loweet therael'inertie of eny unitK 
along the profile (about 1330 TIU), in huge contrast with the adjacent spikes 
which Burk Keyhole Reservoir* 

Line B-B* presents a rather different character in the profiles* The 
Wasatch thermal inertia (fig* 12) is c\ot at ell like that on line A~A'i most 
of its width on line B-B' is an area of unexplainedly low values which mark a 
very distinctive and nonrepresentative area within the widespread formation. 
Wasatch with relatively norma l-appearing thermal-inertia image character 
appears next to the cloud area at the northwest end of the line; there its 
estimated average thermal inertia it 1530 TIU, only 5 percent different from 
that seen on line A-A*. The area of low values does not appear to be related 
to microclimatic factors, nor to any geologic feature of which we are aware. 
For example, neither here nor elsewhere in the image area do thermal-inertia 
values closely and consistently correspond with the inferred llthofacies areas 
delineated in Landsat images* On this line, the Lebo Shale Member of the Fort 
Union has a roughly estimated thermal inertia of 1750 TIU, higher than the 
representative Wasatch values. Most of this is in the broadest part of the 
north-south warm zone, however, so the values almost certainly do not repre- 
sent normal character. The Tullock Member of the Fort Union is estimated at 
1300 TIU, almost 10 percent lower than the Fort Union of line A-A'. Such a 
change is believed to be both real and significant in terms of the geology, 
but no data are available as to possible lithologic changes of the unit be- 
tween the two profile areas. The Pierre Shale has a thermal inertia of 1260 
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TIU, about 5 parcant laaa than obaarvad on lint A-A', and again tht lowttt 
valuaa of any gaologlc unlta on tht prof lit* high and highly varlablt thtnul 
inartla la aaan In tht BlacV Hllla portion of tht profllt at axptcttd in an 
^raa of altarnatlng high vagatatlon danalty and bart rock axpoaurta. 

Fron thla analyalt* It apptara that thtrval data, atptclally whtn coupltd 
with topographic Infomatlon, can aid nattrially In dlacrlnlnatlng gtologlc 
fonatlon and nanbar dlfftrtnctt, aven (at In tht Powdtr Rivtr laaln) where 
unlti art ao varlablt and axpoaurta to poor that gtologiata havt had rtal 
probltna or havt bttn untucctiaful In auch afforta* It apptara that therroal*’ 
inartla dlfftrtnctt of ptrhapi aa llttla at 5 ptrctnt, and ctrtainly 10-15 
ptrctnt, can be dtllntattd and uatd In napping, probably in ttrma of both rock 
units and gtntraliztd aoilt characteristics. So far, patterns seen in ther- 
•al-lnertia Inagts do not natch with vegetation patterns attn in Landsat 
Inages and believed to correspond to subtle facies differences. This problem 
needs further Investigation; it nay relate to difference in resolution of the 
two satellites and to differences in depths ’‘seen” by thermal-inertia measure- 
nents and vegetation root systems. 

2.3. 1.2 Geomorphic domains and linear features 

An interesting and important aspect of using thermal-inertia images is 
that erroneous impressions gained from temperature patterns are corrected and 
a truer picture of surface properties obtained. This is particularly true of 
night thermal data. For example, the very warm sone around the Black Hills in 
the night Inage (fig. 6) disappears in the thermal-inertia image (fig. 7), and 
a whole new pattern energes. The contrast of opposite flanks of the Powder 


37 



tlv«r alfo dlaappaarst indicating no baalc diffartnca in aatarlalt 
acroaa tht rivar* Thus wa can use tharaal, albedo, and tharaal'-intrtla data 
in concert to aaparate physical properties differences (integrated over the 

m 

top deciaeter of the soil or rock profile) fron those effects due to such 
paraneters as slops, altitude, and surface reflectance. 

For the units underlying aost of the Powdex River Rasln, presently avail- 
able theraal data and derived products show aany unexplained patterns, soioe of 
which probably wars due to transient (and now untrackabls) atmospheric 
events. Others, however, are believed to reveal real differences In the 
geologic aaterlals, but current naps In general do not offer a sufficiently 
detailed base for correlation. Certainly some correlations are found with 
Landsat-mapped llthofacles, but more day-night pairs covering varying moisture 
cycles would have been necessary to see through such **noise* as windblown 
sand, surface-wind cooling patterns, and local moisture variations. 

Linear features, often long reaches of streams that appear straight at 
HCMM reaolutlon, are readily defined In the night Images. Many of these 
coincide with breaks or trends In contoured acromagnctic data, suggesting that 
basement tectonic elements have printed through the thick sedimentary sequence 
to control stream coursea. Thla Implies that during aedlmentatlon at earlier 
times, such features effected some control of sedimentary deposltlonal pat- 
terna a conclusion recently elaborated for the Powder River Basin (Slack, 
1981). 

The most remarkable, prevloualy unrecognized, linear feature appears 
prominently on the night Image of 5 and 27 September (figs. 4 and 15) and also 
20 Auguat (fig. 6). Although It la not recognizable as a dlacrete linear 
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fMturc on Landsat inataa (fig* i6a), topographic data (fig* 16b) ahow this 
llnaaaant aa a aubtla drainage divide trending about N55^E* On the thermal 
laagea the a outhward- facing aide ia cooler by 3^ or 4° C» and on the 20 August 
t he raal** inertia image (fig* 7), the aouth aide haa a 17 percent lower thermal 
inertia (1215 TIU) than the north aide (1460 TIU)« Thua, the feature corre" 
latea with a aubdued drainage divide but it cannot appear due to the alope 
effect and muat represent - at leaat in part “ a phyaical property difference 
acroaa the divide. There la no explanation in existing geologic maps (at 
scales f'‘ 0 'i 1:24,000 to 1:500,000) for this feature or why it separates tem- 
perature and topographic domains. The divide is parallel to the prevailing 
wind direction from the west-south-west as shown in eollan deposits south of 
the divide* Moreover, the divide also marks a change in direction of wind 
deposition; deposits to the north are laid down by winds from the north- 
northwest* It is possible that the relatively common eollan sand cover south 
of the divide has controlled drainage habit creating the distinctive topo- 
graphic texture, and the sandy veneer might possibly cause the domain to have 
lower thermal inertia due to lower moisture retention* It is not likely, 
however, that the divide lineament itself is wind related* Extended to the 
northeast, it continues through the linear gap (of the same strike) between 
the Black Hills and the Bear Lodge Mountains* More important, it directly 
overlies one of the most significant breaks (fig* 17) in the aeromagne tic-map 
pattern (U*S* Dept* of Energy, 1979 a, b, and c) of the whole area, it is 
parallel to and roughly coincident with an inflection in the ground-water 
temperatures of the Madison Limestone (fig* 18), and its trend passes through 
several Tertiary Intruslves (fig* 19) and possibly even through Lead, South 
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Figure 16b. — IHumlneted topography Image of the 

Powder River Basin, Nyo., showing the 
drainage divide which is coincident with 
the thermal 1 1 neament . The image was 
computed with a solar declination of 
—7.6 degrees and a local time of 0930 
hrs . 
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Figure 17,— Totel Intenelty megnetio field formllne 
mep, oentrel Powder River Beein, 

Wyo . (U.S.Dept of Energy, 1979e-d). The 
heevy block line ehowe the poeltlon of 
the thermel llnetment. 
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Figure 16. — Ground utter temperttures In tKe Mtdiscn 
Limestone end equlvtlent rocks (Need trd 
others, 1976). Hetvy bltck line Indicates 
position oF the thermal lineament. 
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0«kota location of the HoMstaka Mine. In addition, atreas aMaaureoenta at 
Lead (Aggaon and Hooker, 1960) isply that the preferred direction for norm.il 
faulting would be N30^E, virtually coincident with the orientation of the 
thenaal lineament. Together theae pieces of information provide permissive 
evidence for atructural control relating to this feature. A recently pub- 
lished paper (Slack, 1961) presents substantial corroboration for the con- 
trol. Using subsurface data. Slack proposes that a aeries of northeast- 
trending structural lineaments (fig. 19), which he calls the Belle Fuurche 
Arch, have controlled sedimentation in this area and played an important role 
in determining hydrocarbon accumulation In the southern Powder River Basin. 
One of his lineaments, Cose Butte, Is coincident with part of the thermal 
feature (fig. 19). Additional supportive evidence for structural control Is 
reflected In the shape of the hydrocarbon producing horizons of the upper and 
lower parts of the Muddy Sandstone of Cretaceous age (fig. 19). 

There is further corroborative satellite evidence for our explanation of 
the primary cause of the thermal lineament as a thermal-inertia contrast 
between dissimilar materials. The feature can be seen on two daytime Defense 
Meteorological Satellite Program (DMSP) thermal satellite images (near noon) 
and cannot be seen on a nighttime NOAA-5 thermal satellite image (near 9 
pm). This latter image Is acquired near the time when thermal data should 
largely be Insensitive to thermal-inertia differences because it occurs near 
the crossing times of the diurnal curves. 

We also examined USGS aircraft data across the lineament. Although we 
had not previously recognized the feature on these data, a very subtle thermal 
contrast could be observed in the vicinity of the lineament. Largely because 


46 



of Che elighc ccmperaturt difference end the regional extent of the feature, 
it was not recognieed on the aircraft data, and this reauXt illustrates con-' 
vincingly the potential power of regional thermal satellite data over aircraft 
data for structural'tcctonic analysis. 

A second, parallel lineament, not a stream divide but separating surface- 
textural domains, occurs 30 km to the south (figs. 4 and 13) and also overlies 
an obvious aeromagnetlc break. These features appear to mark fundamental 
structural elements of the southern Powder River Basin and are newly recog- 
nized In HCMM data. To the north, parallel lineaments are marked In the 
images by the Belle Fourche and Little Missouri Rivers. 

Information other than on geologic units per se can be gained from these 
images, most particularly in the demarcation of geomorphic (topographic- 
temperature) domains and in the discrimination of linear features. Night 
images are particularly useful in this regard. The 20 August 1978 image (fig. 
6) provides very clear definition of major areas of distinctive topography, 
commonly linked with distinctive temperature patterns. The north and south 
halves of the basin are distinctly different; the areas east and northeast of 
the Black Hills differ from each other and from the domains of the basin. 
Other finer-scale units also are evident. These do not match the mapped 
geology and, like many geomorphic provinces, are products of a complex devel- 
opment history tied to more factors than the underlying bedrock. We believe 
that important information can be gained in this aspect of the HCMM images, 
especially in understanding surface processes during Tertiary and later 
time. Such information also may lead an understanding of near-surface ground- 
water hydrology across large, diverse areas. We have examined the domains 
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froa tht points of vlow of asppod M^logy, ground*nfst«r ch«Blstry» alntrsl and 
hydrocarbon raaourcaa, and tactonic froMwork, and corraapondancaa are not 
roadlly apparent* This la an aspect of the alaslon data that we did not 
originally anticipate* However, Schneider and others (1979) showed that the 
NOAA (National Oceanic and Aaronautlcs Adolniatration) aatallite VHRR (Very 
High Resolution Radloaeter) data with 900>nB resolution pernltted the defining 
of gcomorphlc donalns very clearly. This work and our present observations 
suggest that careful consideration of satellite thermal Image data by geomor- 
phologists and hydrologists should be undertaken. 

2.3.2 Application to resource studies 

2*3. 2.1 Oil and gas 

It was hypothesized that a few oil and gas fields of the Powder River 
Basin had enough leakage of gas, probably mostly CO 2 , that calcite cementation 
would have occurred near the surface to make the bedrock more resistant, 
possibly to form local topographic highs. A few fields do indeed underlie 
local topographic highs, but this could be fortuitous, and no evidence of 
leakage and cementation has been cited in the literature. However, a recon- 
naissance survey of soil-gas helium shows 37 significant helium anomalies 
within the Wyoming portion of the basin (fig. 20), where most of the oil and 
gas occurs. The reconnaissance scale of helium sampling does not permit an 
accurate comparison with Individual occurrences of oil and gas except for the 
largest fields, but it can be said that all but five of the helium anomalies 
occur over oil and gas fields, and those five are near producing fields. Many 
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oil ond gat flalds do not hava ovarlylng haliun anoaaliaa, and tha inplicatlon 
la conaidarad to ba that aona fialda laak gaaaa upward and nany aora do not. 
Farhapa 15*-20 parcant of oil or gaa fialda night be conaidarad to have leaked 
haliun, although that nunbar night aaally prove twice aa large if nore de> 
tailed aanpling were done. The percentage of fields where gaaea nay have 
produced canentation, aapecially enough to cauae detectable changes in sur- 
face-temperature character, clearly would be expected to be email. Thus, it 
ia encouraging to note that viaual diacrimination auggesta thermal-inertia 
anomalies (20 Auguat image, fig. 7) for 9 out of the 37 areas with anomalous 
helium values (fig. 20). These areas have consistently and noticeably higher 
values than the surrounding areas, in keeping with expectations if cementation 
is locally increased. Most of the areas do not have noticeably different 
topography than their surroundings. Two of the helium-thermal-inertia areas 
do not overlie known oil or gas fields but both are surrounded by fields and 
have several dry holes within the areas. The dry holes may have had only 
subeconomic shows of oil or gas, in which case leakage to the surface voay 
still have occurred; or perhaps holes simply have not been drilled in the 
right places. The two areally largest helium anomalies, one over a giant gas 
field and the other in one over the "barren” areas just described, are marked 
by fairly distinct oval rings in the thermal-inertia pattern. 

Helium anomalies are numerous and areally large in the Montana portion of 
the basin. They occur in a roughly defined ring which corresponds to the 
perineter of a roughly circular area of distinctive topography 73 km wide and 
approximately bisected by the Powder River. This area is virtually without 
oir-gas production and few wells are shown on available source naps. But the 
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ImIIub anoMllts eorrttpond to tht ring of very wira aroot mod In tht night 
20 August laags. Wa know of no raason to axpact such haliua anoaaliat in this 
raglon axcapt In association with oil and gss, and tht wara araas art aostly 
araas of outcrop, parhaps whart caaantatlon Is incraasad> Xnforastlon has not 
yat baan found on facias In the subsurface rocks that alght contain oil and 
gas* If the facies are not truly favorable, uneconomic amounts of oil and gas 
might have been present and leakage could have occurred* If the fades are 
favorable, the area deserves a closer look for exploration* 

An additional point of Interest Is the possible relationship of HCMM 
lineaments and oil-gas occurrence* The major thermal lineaments transecting 
the southern part of the basin (fig* 21) define a block that contains most of 
the significant helium anomalies* Trends in the helium anomalies as contoured 
from present data commonly parallel HCKM lineaments* The two largest fields, 
Fiddler Creek and Clareton, are long and narrow and parallel In trend (and 
between) the basln-trannectlng thermal lineaments* The west end of the Clare- 
ton field ends In a prominent fork, with the southern one of the main linea- 
ments passing through the fork junction* 

2* 3* 2* 2 Uranium 

The uranium districts of the Powder River Basin (fig* 20) have local 
areas of surface alteration as large as 5 by 7 km* These, however, are ex- 
posed dlscontlnuously , and differ only slightly In lithologic character and 
thermal properties relative to the surrounding unaltered ground* If the 
bedrock were totally exposed, an increase In thermal Inertia of perhaps 10 
percent might make the altered ground detectable* In any case, the only 
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th«rnal**lntrtia data aat dota not eovtr tht aaln uranlun araaa; tht 20 August 
data and aast of tht dittrlcta or havt elouda ovtr tht fringt artaa whtrc 
uranium ground occura* Wt havt txamintd low**altitudt aircraft th real data 
for the aitcrtd ground, and to far **nolae" (toll cover, windblown aand, local 
topography, and aolaturt varlationa) attma to completely overwhelm ”algnal" 
related to the discrimination problcm« 

2t3.2.3 Geothermal flux 

Another aspect of our study la to examine the utility of HCMM data for 
geothermal flux mapping* We found that the underground coal flrea now burning 
north of Sheridan, Wyo*, are detectable (barely) on nighttime HCMM thermal 
images. Comparlaon with a mosaic of aircraft thermal images of this area 
(fig. 22) illustrates the scale effects on the appearance of small geothermal 
anomalies. On the aircraft data the anomalies have a sharp, clearly defined 
pattern, whereas the satellite data show an Indistinct pattern which Is not 
distinguishable from geologic and topographic effects. Also, the satellite 
appearance of drainage features is distinctly different from the aircraft 
data. From nighttime aircraft data, the Tongue River appears as a sharp warm 
anomaly with cooler surroundings. The satellite Image, because of Its coarser 
resolution, does not discriminate the narrow water channel, and, thus, the 
drainage appears entirely as a cool zone. 

To examine the expression of geothermal anomalies more fully, HCMM scenes 
of Yellowstone National Park were analyzed. The region has a classic expres- 
sion of moist of the typical hydrothermal features of a vapor^domlnated sys- 
tem. A careful comparison between a nighttime image (fig. 23a) and a detailed 
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Mp of th« hydrothtraol footurot (Saith and Chriotionoon, 1980) (fig. 23b) was 
■ode uolng a boob trantfor acopt. Tht hydrothanal faaturaa coononly ware not 
axpraaaad aa wara anoaalloB (fig* 24a) » and aost of tht vara anomalies on the 
iaagaa art not hydrothtraal-*aatoclattd features (fig. 24b). These results 
demonstrate that nighttime aatelllte Images at this scale are unlikely to be 
useful for detecting similar features elsewhere. Of greater geologic interest 
was the correlation between the caldera outline, some thermal anomaly lows, 
and in particular the anomaly bounded by a sharp edge (fig. 23a) which coin- 
cldes with that part of the gravity field map outlining the southwest side of 
the caldera (fig. 25b). Although this latter anomaly Is a feature with no 
direct counterpart In Landsat Images (fig. 23a), it may provide additional 
Information on the volcanic-tectonic aettlng. From these brief observations 
we conclude that the HCMM data can be useful In understanding the regional 
structural setting of geothermal fields but are not likely to be useful for 
mapping hydrothermal features. 

2. 3. 2. 4 Mapping geologic units In an arid desert environment 

The Cabeza Prieta area In Arizona is an arid desert environment with 
geologic units exposed at a scale suitable for discrimination in HCMM satel- 
lite Images. The area lies very near the major copper district at Ajo and 
contains old prospect developments in hydrothermal ly altered ground. No 
detailed geologic mapping has been previously done of this area as it contains 
a proposed Wilderness site. The main objective of our Investigation has been 
to extend the Interpretation techniques developed In our Powder River Basin 
study. 
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A theraal>intrti« Inagc (fig* 26) was conatructad using tha April 3 and 
4, 1979 scants (AA0342-09150-3; AA0343-20230-1 , 2} with 36-h separation be- 
tween the day and night acquisition tines* The Inage was then compared to the 
AJo geologic nap (Kahle and others, compilers, 197B) and an estlmsic was made 
of the thermal-inertia values of various geologic materials (table 1)* Among 
the sedimentary deposit materials, a wide range of thermal-inertia values was 
found* From highest to lowest values these Included a wet coal-mine dump n«>ar 
AJo (1890 TIU), active pediment areas (1360 TIU), stablilized dunes (1065-1300 
TIU), and active dunes (830-1065 TIU)* Thu most probable explanation for this 
ranking is due to the strong effect of moisture content on thermal inertia* 
Generally, active dunes should have the lowest thermal inertia, as observed, 
owing to their low density and low capacity to retain moisture* 

A somewhat surprlsiiig result was that the thermal inertias of the various 
Igneous rock units were measurably different, indicating a finer discrimina- 
tion capability than previous laboratory data in the literature would suggest 
(Watson, 1979, 1981a)* The literature values of thermal-inertias of Igneous 
rocks show no correlation with either composition or grain size and are indis- 
tinguishable from each other* In Cabeza Prieta, however, we found that the 
felsic Intrusives (together with gneiss and schist) had the highest thermal 
Inertias (>2200 TIU), that extrusive rocks of mafic composition had intermedi- 
ate thermal Inertias (approximately 2000 TIU) and that extrusive rocks of less 
mafic composition had the lowest thermal Inertia (<1900 TIU). We believe that 
the felsic intrusives have the highest thermal inertias because of their high 
quartz content and high surface density and that the differences among extru- 
sive rocks occur because of density differences associated with the amount of 
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Tabic 1*— Estlmatca of thcraal inertia of various geologic aatcrials 


Geologic Material Thermal Inertia (TIU) 

mean range 


Active Dunes 

945 

(830-1065) 

Partially stabiliaed dunes 

1240 


Cluster dunes 

1180 

(1065-1300) 

Active pediment slope 

1360 

(1300-1475) 

Wet mine dump (Ajo) 

1890 


Granite, gneiss, schist 

2200 

(1975-2525) 

Mixed intermediate to mafic 

1750 

(1550-1950) 

volcanic rocks 
Basalt 

1950 

(1700-2200) 
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volatiles prcaent during their foraation and their aaflc content* Although 
these values are suggestive of the possibility of discriainating anong the 
various units, the large overlip in the range of values (hiatograos) indicates 
that classification based solely on thermal inertia differences may not be 
feasible* The overlap between the units can be attributed to several 
causes* These units are exposed in the most rugged terrain of the region and 
thus topographic effects are roost pronounced in these parts of the image. The 
ground that is "sensed" by the technique is a weighted average over the diur- 
nal thermal wavelength (approximately 1 m), with the primary contribution 
coming from the rock and rock fragments of the upper decimeter* Thus differ- 
ences in the effects of weathering processes on the various rock types in this 
environment can produce soroe degree of thermal-inertia differences. 

Because the Powder River Basin was relatively flat terrain the topograph- 
ic effects were not sufficiently important to be considered. In the Cabeza 
Prleta area, however, slopes in excess of 20*^ can be found. To examine the 
contribution of topography, the thermal inertia measured from remote observa- 
tions we considered Che aspects of elevation and slope separately. The eleva- 
tion factor due to changes in the solar and sky radiation, has been examined 
and used to predict the equivalent change in the effective thermal Inertia 
(i*e., that derived from remote rather Chan in situ measurement) as a function 
of elevation (Hummer-Miller, 1981b)* The effective thermal-inertia gradient 
under clear-sky conditions should be roughly 100 TIU per kilometer* In Cabeza 
Prieta, where the maximum relief is 700 m, this maximum correction is only 70 
TIU snd the factor is additive for increasing elevation (i*e*, rocks at higher 
elevations will appear in the image to have lower thermal inertias). 
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The elope factor was evaluated by ualng an algorlthn recently develop<fd 
(Wataon, 1961b)* The turn hai an atlauthal variation proportional to cos(; ” 
37*6^) where la the direction of elope Maaured counterclockwlae from north, 
and thus rldgee with an axial orientation of roughly N40^ will display a 
ttiniroal topographic effect. The noat consistent topographic grain in the 
Cabeza Prieta area is northwest-'southeast and thus generally satisfies this 
constraint* The correction can amount to several hundred TlU's for slopes in 
excess of ten degrees and orientations orthogonal to N40°W (l*e*, N50°E). 

The primary intent of this study area was to examine correlations between 
thermal**inertia values and a variety of common rock types exposed in an arid 
area* In the process of this analysis we also observed, on the night thermal 
image, a number of linear features associated with known major faults (San 
Andreas fault, Garlock fault) and the absence of lineaments in an area near 
Gila Bend which is noted for the absence of structural features (Glia gap). 
We also were able to determine a measure of the satellite's spatial frequency 
response (and thus its ability to detect high-contrast linear features) by the 
observation that lnterstate-10 from Gila Bend to Yuma was observable. 

3.0 MODEL DEVELOPMENT 

Several advances have been made in the development of techniques to 
analyze thermal-infrared data. An algorithm to determine the sensible-heat 
flux from simple field measurements (wind speed, air and ground temperatures) 
has been developed* It provides a direct solution, in parametric form, that 
can be displayed graphically or tabularly. This method has an advantage over 
the previous Iterative solution in that the computation is both very fast and 
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it also provides s cissrtr understanding of the drag coefficienti with its 
variation and response to different conditions. At low wind speeds the drag 
coefficient cannot be treated as a constant. Both the conputationsl speed and 
analysis of the drag coefficient can be importa.it for rcmotc-’scnsing applica- 
tions involving thermal scanner data (Watson, 1980). 

A substantial advance was the development of a method, based solely on 
remote-sensing data, to estimate those meteorological effects which mubt be 
knovm for thermal-inertia mapping. It assumes that the atmospheric fluxes aru 
spatially invariant and that the solar, sky, and sensible heat fluxes can be 
approximated by a simple mathematical form. Coefficients are determined from 
a least-squares method by fitting observational data to our thermal model. A 
comparison between field measurements and the model-derived flux shows that 
good agreement can be achieved. An analysis of the limitations of the method 
was also made (Watson and Hummer-Miller, 1981a). 

This new method of estimating atmospheric parameters was the basis for a 
revised thermal-inertia algorithm (Watson, 1981a). The new form is: 

- ("p . ^ V + C(* ,i) (A-A^j))A‘.Vij 

where Aj^j and are the corresponding albedo and temperature difference of 
the ith pixel and jth line. A and V are the mean values for the area in 
question, and P is a select value for the mean thermal inertia (generally I5uu 
TIU). C(x , 6) Is a function of the site latitude, X, and the solar declina- 
tion, 6. The advantage of this algorithm lies in the fact that we are dealing 
with albedo and thermal differences rather than absolute values. Thus, the 
computed thermal inertia Is less sensitive to offsets caused by calibration 
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•rrori or ataotphcrlc backoeottorlng and cranaalaaion affacta* 

Othar ■odtling atudlai cantorad around davaloplng an algorltha for alava** 
Cion corractlon of taapcratura and tharmal~lnartia inagas (Huat&ar-MlUer, 
1981b). Thay ara baiad on application of tha linaarlzad Fouriar aarlat method 
(Wataon, 1975; Wataon» 1979) to almple forma of tha aolar flux (Hummer-Miller » 
1981a) derived from a rapraaantatlva aat of field obaarvationa. It was found 
that flux varlatlona with elevation can cauae changea In tha mean diurnal 
tamparatura gradient from -4° to **l<i° C par km (evaluated at 2000 m). Changes 
In the tamparatura-dlffaranca gradient of par km ara alao produced and 
these are aquivslent to an affective thermal-inertia gradient of 100 TIL) per 
km. 

In addition, a simple topographic slope corractlon method has been devel- 
oped using the linearized thermal model and assuming slopes leas than about 
20°. The corractlon can be used to analyze individual thermal images or 
composite products such as temperature difference or thermal Inertia. Simple 
curves were determined for latltudeb of 30° and 50° (Watson, 1981b). The form 
Is easily adapted for analysis of HCMM Images using the DMA (Defense Mapping 
Agency) digital terrain data (Watson, 1981b). 

A major concern In this Investigation has been the accurate registration 
of day and night images. We have developed an Image-registration algorithm 
which appears to be substantially better than the current registration prod- 
ucts provided by NASA (Watson and others, 1981). The initial test of this 
algorithm used the 20 August 1978 data of the Powder River Basin. A small 
number, less than ten, of very clearly delineated features, generally the 
water-dam Interfaces of reservoirs, were selected as control. Subsequently, 
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«n affint trans fora* cion waa dataralnad by bait fitting thaaa points* Our 
first taat Indicataa a rtaidual arror of < 2 plxala. Iha NASA product for thu 
aana acana dlaplaya arrors of aany plxala rasulting in "doubla drainage" 
*affacts and an offaat of aavaral kiloaatara. 

During tha initial atagas of our axparlmantation with control points, we 
axpactad that the drainage pattern in the Powder River Basin was substantial 
enough to pro> da extensive control for registration. We discovered, however, 
that drainages are often unreliable identification features and the resulting 
control points were too Inaccurate to provide the transforoation coordi- 
nates. We also experimented with cross-correlation techniques in the Cabeza 
Prieta area but were unsuccessful owing to the strong topographic grain. The 
affine traneformar ion which we employed has the additional advantage that it 
can be adapted to very fast computer processing schemes (Braccini and Marino, 
1980). 

4.0 DIGITAL IMAGE PROCESSING 

This section presents an outline of the image processing techniques used 
in this study. Figure 27 is a simplified flow diagram of the basic processing 
steps; the computer programs referenced in this diagram are Included in Ap- 
pendix A. The initial processing Involves obtaining the HCMM computer- 
compatible tapes (CCT) and altering the data format to be consistent with our 
computer software* Sometimes the area of interest spans two scenes and, 
consequently, must be appended into a single file. To produce enhanced images 
of these products, the appropriate area of the image is statistically sampled 
to form a histogram and to derive the mean, median, mode, variance, standard 
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dtvlation, and cumuUtlvt fiaquanclts. With thaaa •tatlstlci a decision Is 
■ade as to how the scene contrast ahould be enhanced using an appropriate 
transformation of the scene brightness. From our experience a 1 or at most 2 
percent linear stretch of the data produces a good starting point* Thus the 1 
or 2 percent and 98 or 99 percent points are transformed linearly to the 
extremes of a 256-step gray scale (0 and 255) and the center density value 
(Dn) of the distribution Is transformed linearly to the center of the gray 
scale (127 Dn). Other useful products Include linear stretches on paper where 
a ingle print character Is assigned for each Dn value and color-coding on 
f lxk3* 

The next processing operation is performed to register the night-therma] 
file to the day files. This operation, including a oiscusslon of the general 
considerations, are detailed in Watson and others (1981). Control features 
are selected from the positive transparency images produced from the CCT and 
locations arc measured to one pixel accuracy. These values are then used to 
determine the affine coefficients for a best-fit transformation (REGALG). 
This transformation provides a rotation correction for the inclined satellite 
orbital tracks, an origin shift, and scale changes both along and across the 
scan line. The actual registration of the night file is then performed using 
the affine transformation coefficients in the GE0MX4 computer program. This 
program assigns radiometric values to the newly registered image employing a 
nearest-neighbor method. 

At this stage the data are In the appropriate format for thermal-^inertia 
and temperature-difference mapping. The algorithms which we use (Watson, 
1981a) employ average scene values of albedo, day temperature, and night 
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t«mptrature> The portion of the scene from which these vslues sre determined 
Is bssed on the sssuoptlon of stmospherlc Invsrlsnce end thus, ss s minimum 
eonstrslnt must be cloud-free on both Images. The program RAPFIT is then used 
to compute the appropriate coefficients for this algorithm, and the program 
HCMTIDT is employed to construct both temperature-difference and thermal 
Inertia files. The resulting files are then processed using these techniques 
described In the beginning paragraph of 4.0. 

Another analysis technique used in this study Is to construct profiles 
across the Image data. The end points of particular profiles are chosen and 
digital values along the line are obtained for all products: day reflectance, 

day thermal, registered night thermal, and thermal inertia using a nearest- 
nei^rbor algorithm. The corresponding elevation profile Is obtained by digi- 
tizing the appropriate portion of a 1:250,000 topographic map and adjusting It 
to match the satellite data. This task is made easier if the profiles cross 
distinct features such as reservoirs and rivers. After the elevation data are 
registered to the satellite data, the p' ogram PROFLE is used to plot the 
profiles and cross plot pairs of data values. The profiles can be plotted 
various scales and thus directly overlaid on any base material for compari- 
son. The cross-plotting option is valuable for examining correlations (for 
example, we observed that the day thermal versus elevation data fit the adi? 
batlc lapse rate). 
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5. 1 APPENDlXEii 


l^iiiting* of tht following computer programe described in section 4.0 
follow: 


REGALG 

RAPFIT 

GEOhCU 

HCMTIDT 

PROFLE 
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10 I mCALC To comf^uit te»t fit «fline tran»lorir*tion for HCMH rtgfst. 

30 I rrogrtm language is eitenoca Basic. 

10 I Biogism written by K.nsison 

40 ' OPTION BASi: 1 
00 STANOARD 

*0 INTEGER 1 .Pr inter 

70 DIM Aoiu(lOU) .yoloOOO) .Anew(lOO) .ynewnoo) .Nsmc^lOO) , S ( 3) , R ( 3) , XI ( 1 00) ,y t O DC) , 1 fe: ) , C(. jn 

(}U0) ,PmI <4) , E«) (4) 

BO PRINTER IS lb 

BO LlNPUr *Oo you wish to r eso sp input file y pi enter vslucs ir.snusHy N*.Anz5 
100 If Ani$«*y' THEN GOTO 130 

110 LINPUT *Do you wish to store an input tile y or N'.AnyS 
120 IF AnyS-’N* THEN GOTO 350 

130 IL1NPU7 "Enter Filename lor Day(new) point s* ,r i lelS 
140 IF Any5«'y THEN CREATE FilelS.lO 

150 LINPUT "Enter Filename lor Nite(old) points' ,Fi le2$ 

HU IF Any5«"V" THEN CREATE File2S,10 

17U IF Ar.y$«*y THEN GOTO 350 

IdO ASSIGN |1 TO FilelS 

190 READ (1;N 

200 Nnum-N 

210 REDIM Counter <N) 

220 READ |1;G$ I Implies a title which is ignored 

230 FOR 1-1 TO N 

240 READ «} iVnewd ) ,xnew(l) 

250 Counter(I)-! 

260 NEXT I 

270 ASSIGN tl TO FilelS 

260 ASSIGN fl TO File26 

29u READ tl;Jmax 

300 READ flrCS I Ignored title 

310 FOB I-’l TO N 

320 READ il ;yolo(l ) .Xold (1 ) 

330 NEXT I 

340 GOTO 530 

350 INPUT "Enter numoer o£ control point quadrupl es" ,N 

360 REDIM Xolo(N) ,yold(N) ,Xnew(N) ,ynew(N) ,Xt(N) ,yt(N) .Counter (N) 

370 Printet-16 

3d0 PRINTER IS Printer i For debugging 

390 LINPUT "Enter a descriptive name <50 char lor this r<»g i strat ion* , Name $ 

400 INPUT "Enter tne maximum lines and pixels trom old im.ag e* .Lrnax , Pmax 
410 PRINT PAGE 

420 PRINT LIN(IO) ,SPA(5) ,*To enter control points use the formats* 

430 PRINT ■ Day(new) Night(old)" 

440 PRINT * scan no, pixel no , scan no, pixel no* 

450 FOR 1-1 TO N 

460 I NCTCsScan-Line 

470 PRINT L1N(9),SPA(20), "Control po i nt : * ; VALS (1 ) 

460 INPUT "Scan new ,P i xel_new, Scan_old , P ixel_^old" , yncw(l ) , xnew ( I ) ,yol d{ I ) , Xol d ( I ) 

490 Counterll)»T - - ~ 

500 PRINT PAGE 

SIO NEXT 1 

520 Enter: PRINTER IS It- 

530 PRINT PAGE 

540 PRINT LIN(IO) ,SPA(16) DATA CHECK* 

550 PRINT "Control no NE*i OLD" 

560 PRINT • Line Pixel Line Pixel* 

570 FOR 1-1 TC N 

5B0 PRINT USING 590; Counter ( 1 ), y new ( 1 ), xnew) 1 )» Volo ( D «xo Id 1 1 ) 


SVU IMAGE POD, 6X,ODOODOO. 4X,DDOOODD« 6Xr DODDDDD .6X. OODODDO 

fcoo 1 

ftIO LlNPUT *Af« values coitect -y or N 7”,Ans$ 

§20 ir AnsM)>l)**y* THEN GOTO Cont 

§30. INPUT *Do you WISH to otlctc values ycs or NO 7" ,Acj5 

§40 IF Aq$ll,l|0*y* THEN GOTO S4U 

§b0 INPUT *Hoh many controls to delete 7*fNurr.oer 

§60 DIM Check(SO) 

670 REDIK CnecK ( Numoer ) 

660 PAINT 'Enter control numbers in any order* 

690 MAT INPUT Check 

700 KiAT SORT ChecAC) EES 

710 FOR 2"1 TO liumoer 

72u N»N-1 

730 N^T SEARCH Counter (•), bOC (CnecA ( Z) ) 1 1 Star t 
740 FOR l«IStat t TO N 

760 Xololl) >kola({4l) 

760 yolo(l) >yolo(l4l) 

770 Xnew(l) ■Xnew(l4l) 

780 yni>( 1) «ynew(l4l ) 

790 Counter ( I ) -Counter (l + l ) 

600 NEXT 1 

610 NEXT 2 

620 RXDIM xold(N) ,Xnev(N) ,yold(N),ynew(N) 

630 GOTO 530 

640 INPUT 'Enter control number for incorrect values', K 

650 PRINT LIN(20) ,SPA( 5) , ' V new , xnew, V ol d , Xold- * j V new( K ) ,xnew (X) , Void (K ) , Xol d (K) 

860 INPUT 'Enter correct values lor all 4 par atr,e ter s* ,ynew ( X) , Xnew (x) , y ol d(K ) , Xold (K ) 
670 GOTO 530 

660 Cont; LINPUT *Do you wxsh to list control points on pr i nte r .*, Lists 


890 

IF LlstSU ,l)0'y* 

THEN 

GOTO 

960 


900 

PRINTER IS 0 





910 

PRINT 'Control no 


NEW 


OLD' 

520 

PRINT * 

Line 


Pixel 

Line Pixel* 

930 

FOR 1-1 TO N 






940 PRINT USING 590 ; Counte r (I ) ,y new ( I ) , Xnew ( 1 ) , yold( 1 ) , Xo Id ( 1 ) 

950 IMAGE DDD,7X,DDDDD,2X,DDDDD, 9X, DDDDD , 2X , DDDDD 

960 NEXT 1 

970 PRINT L1N(5) 

S60 CALL Bxlinear(Xnew(*) ,ynew(«) ,xold(*) ,R(») ,N) 

990 CALL Bilinear{Xnew(*) ,ynew(*) ,yold(») ,SC) ,N) 
lOOC I wnere: Xold«R(l ) *xnew4R (2) • ynew+R( 3) 

1010 1 yold«S{l)*Xnew+S(2)*ynew4S(3) 

1020 1 Thus to transform a niqht image the new DN values are computed: 

1030 I DNnew(PiAel,Line)-DNold (P,L) 

1040 I where P»R( 1 ) *P ixe 1 +R ( 2) 'Ll ne+ R ( 3) 

1050 I and L-S (l)*PiAel4S(2)*Line4S (3) 

1060 PRINTER IS 0 

1070 Dell2»R(l )'S(2)-S(1)*R(2) 

1080 Dell3-R(3) 'S (1) -S (3) ‘Rd ) 

1090 Del23-R(2) *S(3)-R(3)*S(2) 

1100 Pxl (1)«S(2)-R(2) 

1110 Pxl (2) >S(2)*Prr.ax-R(2) 

1120 Pxl(3)»5{2)'Pmax-R(2) »Lmax 

1130 Pxl (4) -S(2)-R(2) *Lmax 
1140 Lxl(l)-R(l)-SU) 

1150 Lxl(2)-R(l)-S(l)*Prr.ax 
1160 Lxl (3)-R(l )»Lmax-S(l)'Pmax 
1170 Lxl(4)-R(l)'Lmax-S(l) 

IIBO MAT Lxl -Lxl+lDell 3) 


IS 


ORIGIMAI PAGE \S 
OF POOR QUALITY 



mo 

1200 

1210 

1220 

1230 

1240 

12S0 

1260 

1370 

13B0 

1290 

1300 

1310 

1330 

1330 

1340 

mo 

1360 
1370 
1380 
1390 
1400 
1410 
1420 
14 30 
1440 
14 50 
1460 
1470 
1460 
1490 
1500 
1510 
1520 
1530 
1540 
1550 
1560 
1570 
1560 
1590 
1 600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1660 
1690 
1700 
1710 
1720 
1730 
1740 
1750 
1760 
1770 
17B0 


HAT Ul>Lil/(Dcll3) 

MAT Pal«P«l« (Ocl23) 

HAT Pil>P«l/(Dcll2) 
hat search Pal (•) ,MAX{ Parr.aa 
MAT SEARCH Pa 1 ( • ) ,H1 N i Pa ir.t n 
MAT SEARCH La 1 (• } ,MAX; Lamaa 
MAT SEARCH La 1 (• ) , Ml N ; Lamn 
Dpia«Pair.aa>Pamtn4l 
Dl in*Lairaa - Lamin4 1 

Alph**R(l )• ( Patr,in-1 )4R(2) • (Lamin-1 ) *R( 3) 

Be l»»S(l )• ( Pamin-1 ) ♦$ (2 > • (Lamin-1 )4S(3) 

C2*Alpha*S(2)~Bet«*R(2) 

C3*Alpna*S(l)-Bcla*B(l) 

Pia_old_inean«SUM( Xold) /N 
P ia”new”mean»SUM( Xnewj /N 
Lin”new”riean»SUM ( y new) /N 
Lin'“olo”niean"SUMiyold)/N 
Drl>Dell2 

Pia_oil»s (2)A>el*Plx old_mean-R (2)/Del*Lin old_tnean-C2 /Del-Pia new_mean 
Lin”oil»> (S ( 1 )/Del ) "Fia old iriean+R (1 j/Del'Ixn old mean^C3/De 1-Cin__new_rican 
CoE”tneia»(R(l)'*S{2) )/2" 

Sin“tnetaw(S{l)-R(2) )/2 

deg" 

Theta»ATN(Sin theta/Cos theta) 

Mag k-SwR(R( 1T'2+S{1)‘27 
Mag"y.SOR(R(2)*2+S(2) *2) 

PRINT •SOLUTION FOR!* 

PRINT Names 
PRINT 

PRINT "Original image is •jPrr.ax;* pixels by “;Lir.ax;* lines. * 

PRINT 

PRINT 

STANDARD 

PRINT •R»";R(1) ;R(2) ;Alpna 
PRINT 

PRINT *S»*;S(1)!S(2) jBeta 
PRINT 

PRINT ’EK (toc:R(3)» *!R(3)j* S ( 3) • ’jSO) 

PRINT 

PRINT *Rhere: DNnew (P ixe 1 , Line) -DNold ( P , L) * 

PRINT 

PRINT * P«R( 1 ) •Pixel4R(2)*Line*R(3) * 

PRINT * L»S(1 ) •Pixcl+S(2) •Linc4S(3) • 

PRINT 

PRINT * New image size: *;Dpix;* pixels by ■;Dlin;* lines.* 

PRINT 

PRINT * Mean rotation angle :*;Theta 

PRINT * Magnification s-X *;Mag xj* -y ";Mag__y 

PRINT 

PRINT * Offset! Pixels *;Pix off 

PRINT * Lines *}Lin"off 

PRINTER IS printer 

Plotting: CALL Plo t {* I ni t ial contr ols: * , Xold (• ) , yold ( • > , Xncw( • ) , ynew( • ) , n , Nane5 ,C c^r. tc i l * ) 
Residual-0 
PRINTER IS 0 
Max_rcsidual-0 

PRINT * Control no. Yt Xt Residual 

Del-Dell2 

C2- Del23 


17f0 

moo 

mio 

1120 

1S30. 

2«40 

1I»0 

moo 

U70 

moo 

moo 

1900 

1910 

1920 

1930 

1940 

19S0 

1 960 
1970 
I960 
1990 
2000 
2010 
2020 
20 30 
2040 
2060 
2060 
2070 
2060 
2090 
2100 
2110 
2120 
2130 
2140 
2160 
2160 
2170 
2160 
2190 
2200 
2210 
2220 
2230 
224 0 
2260 
2260 
2270 
2260 
2290 
2300 
2310 
2320 

2 330 
2340 
2360 
2360 
2370 
2360 


C3>Dtll3 
roli 1*1 TO N 

*t<l)-(Xold(I)*$ l2)-yold(n*R(2)-C2)/Del 
ytm*~(Xold(l)«S(l)-yold(n*IMl l-cil/oel 
Han I fd-txHl )-Xnew(l) )*24(yt(n-ynew(n ) *2 
FMNT Counlcc (1) . Xt (n ,y t(I) , SUR (Mai_,c cd) 

X«aidual*MaR_rtd4Rraidual 

MaR_:aa>oual*KAX (Hai_rcd,MaR_iei tdual ) 

ir Max (csidual*Max Fto TriCN"lR>ax rts>Counter (1 ) 

NEXT l" " " 

Nax_^rcs>SOR(Max_rca>oual 1 
RcaToualaSJRiReaioual/iN-l )) 

FIXED 2 
PRINTER IS 0 
PRINT 
PRINT 
FIXED 2 

PRINT * Mean icaioual vector lcn 9 th ■* {Residual;” pixels* 

PRINT * Max residual vector length ■*;Max res;* pixels*;” at point * ; 1 rrax res 
PRINT 

CALL Plot (*T tans formed controls; ”,Xt(*),Vt(*) ,xnew(* ) , ynew(* ) ,N, Name 5, Counter (• ) ) 
GOTO Enter « 

END 

SUB Bilinear (X(*) ,Y ( *) , 2 (* ) ,M {• ) ,N) 

I Z»M(1 ) *X4M(2) *y4K(3) LEAST SQUARES ESTIMATE FOR M 
OPTION BASE 1 

DIM S(100,3) ,Transpose(3J00),Inverse(100,lC0),DuiT(100,100) ,Zl(100tl) ,M1(3>1) 
REDIM S(N,3) .Transpose (3 ,N) , 1 nvei se { K , N) ,Dum(K,N) ,23 (N,l) 

FOR 1>1 TO N 
S(l,l)-X(l) 

S(i,2)»yu) 

S(1 ,3) -1 
21(1,1)»2(1) 

NEXT 1 

MAT Transpose»TRN(S) 

MAT Dum, “Transpose'S 
MAT Inverse*INV(Dum) 

MAT DumBl nver se*T ranspose 
MAT Ml “D urn' 21 
FOR J-l TO 3 
H(J) -Ml (J,l) 

NEXT J 
SUBEND 

SUB Plot {Plot title?, Xnew (•) ,ynew(* ) , Xold (•) , yold(* ) ,N, Name?, Counter (• ) ) 

OPTION BASE 1" 

INTEGER I 
STANDARD 

PLOTTER IS 13 ,*GRAPH1CS* 

GRAPHICS 

LOCATE 5,160,6,95 
MOVE 10,97 

LABEL Plot titleSs* *SName? 

MAT SEARCH"Xnew,MAX;Xmax 
MAT SEARCH Xold ,MAX; Xmaxl 
Xroax-MAX (Xmax , xmaxl ) 

MAT SEARCH Yol d ,MAX; Ymax 1 
MAT SEARCH y new, MAX ; ymax 
ymax -MAX ( ymax , ymax 1) 

MAT SEARCH Y new, MIN ; ymi n 
MAT SEARCH Void, MIN; Yminl 
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2310 

2400 

2410 

2420 

2430 

2440 

2450 

2460 

2470 

2460 

2490 

2*)00 

2S10 

2S20 

2330 

2340 

2330 

2360 

2370 

3360 

2390 

2600 

3610 

2620 

2630 

2640 

3630 

2660 

2670 

3660 

2690 

2700 

2710 

2720 


ynin*MiN(ymln 
HAT SEARCH Xo Id ,M1 N t Xmt nl 
MAT SEARCH Xn*w ,M1N< Xmi n 
X»in-MlN(xmin,Xinin) ) 

SHOW Xmln- 33 |Xmax« 23 ,ymaK 433 |y rin -23 

UNE TyPE 1 

FRAME 

LORG 2 

CSIZE 3. 3 

FOR 1-1 TO N 

MOVE Xo)d(I ) ,yold(I ) 

DRAW Xncw(l ) ,ynew(I ) 

Xx-Xnew(l)-Xold(l) 

yyynew(X)-yolod) 

IF yy «0 THEN GOTO 2330 
LDIR Xx.yy 

MOVE xnew(i ) ,yn«w(X) 

XF Xx >«0 THEN LABEL *>* 6 VAL 5 (Counter ( 1 ) ) 

XF Xx >>0 THEN GOTO 2620 
LORG 8 

LDIR -XX, -yy 

LABEL VALS(Counter(X) ) 6 *<* 

LORC 2 
NEXT I 
SETCU 
LDXR 0 
CSIZE 3.3 
MOVE 69,3 
LORC 3 

LABEL ’COST to continue! DUMP GRAPHICG for h»rd copy* 
PAUSE 

EXIT GRAPHICS 

SUBEND 

END 




1 ^ 



10 I OAFriT To ptitora tht itoloncl ctmoiiphtiic (ittino, 

30 I rioorcm lanoutot is tuttnOcd 0*tis. 

30 I ffooiam wittttn by K.WAtton 

40 I Omapptng Algor tthm dtvtlopmrnt . 

30 i ENTER; SITE LATITUDE AND SOLAR DECLINATION. 

40 I RROCRAM COMPUTES Td«y AND Tnttt FOR THE HCPU<t TIMES , 

70 I VARY INC THE TtAy AND Sky t«Cior. 

SO I INPUT OBSERVATIONAL DATA? 

• 0 I Tdayo.Tntito; Albroo 

100 I COMPUTE THE St«nd«cd Devi«tton for te^p d«t« - SriAde print output. 
1)0 Piconi>SwR(PI*34C0*24) 

130 OPTION BASE 1 

130 DIM rluKl ( 24) ,Temp<34) .Phi (34 I .rlui (34) ,S tt lr«i S( 10) I ftO] 

140 DIM Pl(34 ,24) ,P<>(24.24 ).Cosx(24) ,Dcv(20 ,20 ) 

140 INTECEK 1 ,K,L,Oel Ik 
UU DEC 

no I RECALL THE Phi . 

160 ASSIGN |1 TO ’Phi 

190 READ IlsM 

200 IP K>24 THEN OOTO 330 

210 PRINT ’ERROR I j’jMj’ P VALUES - SHOULD BE 24’ 

220 STOP 

230 READ Hi Fhi(*) 

240 ASSIGN |1 TO ’Phi_m’ 

340 INPUT ’Enter the kite la 1 1 tude ( Deg tees) *, Latitude 
260 INPUT ’Enter tne aolat decl i na t ion (Deg re ec ) ’, Solar aec 
270 OUTPUT 9j’R’ 

260 tl.TER 9jTS 

290 R-FI/12 

300 wnFl/(12’360Q) 

310 SO-1360 

320 Sigra«4. 67E-6 

330 I FIXED PARAMETERS. 

340 Eric&ivity«l 
350 Slcpc-0 
360 Azinutri-O 
370 Flu»_oi£ r ute»0 

36u INPUT ’Eniet tne mean aay Tetrp (DEG K) iiom scene data’ ,Tdayo 

39U li.PUT "Enter tne mean nite Temp iDEG K) fiotr scene oata’.Tniteo 

4 00 Cl -CCS (La ti tode) "COS (Solar_oec) 

4l0 Si *£1 N (Lat i tuoe )»S1N (Solar ^oec) 

420 CO£Z ( 1) -Cl+Sl ” 

4 30 For. K-2 TO 24 

440 Cesz (K) -Cl'COS ( (R-l)*15)4Sl 
450 II Ccsz(K)<0 THEN Cosz(K)-0 
460 NEXT K 

470 CzJ-SlN(Sl&pc) •SlN(AzimutH) *COS(Solar dec) 

46 0 CZ2-COS (La ti tuoe) ’COS (S lope ) -SI N (S lope ) ’COS (A z i mJ t h) • SI N (La 1 1 tude) 
490 Cz2 -Czi’COS (Solar dec) 

40 0 Cz3-blN (La ti tuoe)'’COS (Slope )+COS (Latitude )»CCS (Az in-uth) -SIN (Slope ) 
410 Cz3-Cz3*S JN(Sclar__oec) 

420 INPUT ’Enter mean^A.' .Albedo 

430 Inertia-1400 

440 Oeviatior.2-1 0000 

44U Ccns2 -1 net t ia/T icon* 

460 I LOOP parameters. FIRST LOOP 

470 PRINTER IS 0 

460 PRINT • Tsky’.’Sky factor", * Deviation* 

490 Return: PRINTER IS TS 
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*00 lUFUT 

*10 INPOT "T^fcy'ptn" ,T»«y“iMn 

*20 INPUT *Oit«y* »Dt*ky 

*20 INPUT *i«y (pcioi m«ii*(|ky f«ctoi min 

*4U INPUT •Sky"l *cior"»»n* ,iky”< *cloi“m»n 

**0 INPUT *0»ky 1 «cto7* iDiKy^i Zciot 

**0 Dtv minalCl'* 

*70 Ktni Tor imy-TkKy^iRtn TO T»liy^»*« STEP Dtsky 
*10 PRINT •Tiky»*|T*«y 
*t0 riuR SRy«Si9mi*Ttky *4 

700 rOR 5ky l«ctor>SRy t«ctof mtn TO Sky (actot max STEP DBRy factor 
710 OUTPUT 5|*R" " " 

720 ENTER tjMonth.Oay, Hour , Mi nutt .Second 

720 PRINT TAS(20) {"Sky factoia*|Sky factor, Hour;* t *;Minute 
740 riu« ditkCt«So* ( I'Sky factor) 

750 rOR T»1 TO 24 
760 X«(l-l»*15 

770 Coszl ■COiX(l ) 

7S0 ir CoaiKO THEN Cotil*0 

790 rluxl(l)a{l-AlD»oo)*(riUk Oif fuae^Coazl'Flux dircct)4i:r.iBiivity*rlux aky 
*00 NEXT 1 ” 

*10 SuraO 

*20 FOR 1-1 TO 24 

130 Sur-Fluxl (l)4Surr. 

440 NEXT 1 
450 Suk-Sun/M 

460 Vo- (5uiT/(EiriBaivity*Si^:<''>a))*.25 

*70 Alpr.a-2*Err.iB8lvity*S s«^(r,a*Vo‘3 

*a0 Conx-Ccnx2/A Ipna 

490 Cons) -Aipha*4/3 

900 FOR 1-1 TO 24 

910 Flox(I) -Fluxl(l)*Alpha*Vo 

920 NEXT 1 

930 FOR X-1 TO 24 

940 FOR K-1 TO 24 

950 L-I-N-rl 

960 IF Ul THEN L-U24 

970 Dcllk«0 

960 IF I-K THEN Delix-1 

990 P5 (1 , K) -Cor.t2*Pni (L) -f Cons) ‘Del ik 

1000 NEXT K 

lulu NEXT 1 

1020 MAT P1-1NV(P5) 

10 30 J FOR.'I Ter.p 
1040 MAT Teir.p-Pl* Flux 

1050 I REMEHdER Te/rp(K) is at ti.ua K-1 nour s. {r.casurad frotr. noon) 

1060 Toay(Tc.7,p(2) 4Tair.p(3) )».5 
10 70 Tni te- (Tamp (14)4Terrp(15))*.5 
lOoO 11 - (Tsky-Tsxy_^min)/Dtsxy4l 

1090 Jl»(Sky factoF-Sxy i actor min)/Dsky_tactor4l 
1100 Oav( 11 , j1) -5uR( (T cay-Tday o) * 2+ (Tni ta-Tni teo) * 2) 

1110 Dav nin«MlN (Dev_r.in,Dev (11 ,J1 ) ) 

1120 IF Dav ( 11 ,0 1) -Dev min THEN Tsky opt-Tsxy 

1130 IF Dev (1 1 ,J 1 ) -Dav"min THEN Sky Factor opt-Sky factor 

1140 PRINTER IS 0 " 

1150 PRINT Tsxy.Sxy f actor .Dev (1 1 ,J 1 ); 11 ;J1 
11*0 PRINTER IS 16 ~ 

1170 NEXT Sky factor 
11*0 NEXT Tsky 
1190 BEEP 




OMQINM. PMX It 

or MXMqUMJEt 


Uwu PM1N7 *OpttmuRi lili Tft»y •*jTt«y opt«*l>ky (acioi ■*)S«y itctor opt 
IdiO GOTO Pttuin " • “ ” ” 

1220 INPl/ 1 ‘intti THE 1 (Tft«y} J (S^y iDCtor) i/ALUCS”« 1,J 
1230 Tfc*yo»3J0* (lo-l)»lO 
1240 S«y I tcturoa ( Jo-1 } * .U& 

12‘!>J I “ LOOP PARAKETtPS. SECOND LOOP 

1240 FOR Tl.k>-TtftyO-2 TO T»Kyc«2 STEP 1 
1270 PRINT ‘‘7tK)«*{TftNy 
12bi) riuk kk)-«S > jn «*Tt«y * 4 

12»o fvh 5».y i*ciui*SAy i«ctorO'.02 TO SRy iActoM.02 STEP .0) 

J3oO OU'IPuT ?(*A' 

1310 Ei>7i.P d!Mcntn,Oay .Hour .Ktnote ,St'Cond 

13^0 PR1..T TAfc(20) {’S»y_l ACiof •";S*iy ( *c lor , Hour i * i •• j M 1 n uie 

1330 Flu* Cl 1 «ct«Sc* .k»"j -S»y lactotT 

1 340 FOR T»1 TO 24 

13SD XMl-l)»lb 

1 3t0 Cctzl aCosz ( 1 ) 

1 370 IF CotzXO THEN Ccttl»0 

1380 Flu«l (l)>(l-Alotoo) * (Flux if 1 otc«Cotzl *Fl ox d ir ret) 4Emi cti v i ty * Fl ox txy 
1380 neat 1 ” “ 

1400 SoPaO 
1410 FOR ]a) TO 24 
1420 Sur«riuxl (I ) 4SUIT 1 
1430 NEXT 1 
14 40 Sura£ur/M 

) a£>o VC- iSur,./ (Lrimvity'S igrra) )* ,2S 
I4ti0 Air'r>a*3*E;Mxxivity*5 igma* Vo* 3 
1470 Con&-Cor.c2/Alpn* 

I4b0 Cont)aAlpr.«*4/3 

1480 FOR lal TO 24 

ISoO FI ux ( 1 ) *ri ux 1 ( 1 ) 4A Ipna* Vo 

1610 NEXT 1 

1620 FOR lal TO 24 

1630 FOR K-1 ro 24 

1640 Lal-K4l 

1660 IF L<1 THEN LaL424 

1660 OeliNaO 

1670 IF XK then D*1iX-1 

1660 P6|l,K)“Ccn62*PM(L)*Consl»Dellk 

1690 NEXT K 

1600 NEXT 1 

1610 MAT PlalHV(P6) 

1620 J FORK Ter.p 
1630 MAT Teff.papl -Flux 

1640 1 REMEMBER Tr(ip(K) It At tlire K-1 hoort. (rt.cASured iron, noon) 

U6U Tcay- (Terp{ 2) 4Tcnp( 3) )• .6 
1 o6o Tni le »(Tcrp (14) 4Ten.p( 1 6) ) • . 5 
1670 ll-(Tsky-260)/54l 

1660 J1 »Sky_i ic tor/. 0 6 + 1 

1690 Dev ( IlTo 1 ) »5wR( (7oAy-Toay o) ‘24 (Tnite-Tni teo) *2) 

1700 PRINT Ttxy.Sky f AC tor ,Oev (11 . J1 ) 

1710 NEXT Sky 1 ACtoT 

1720 NEXT Ttxy 

1730 Min devajcis 

1740 FOR'la) TO 5 

17 60 FOR J a) TO 6 

1760 Min oevar,lN(Min dev , Dev ( l , J ) ) 

1770 IF Dev (1 ,J) aKin"dcv THEN loal 
1760 NEXT J 



179U NlXT 1 

JsOO T»^y« X»Xyo*( 1 0-3) • 1 

3fcli> Siy ^«ciot*SKy t«ctoro4( jo-3)* .01 

ls2U STOP 

lo30 CHO 


n 



PR0GRM1 0E0nX4 


t 

C 

C RCnOTE SENSING ARRAY PROCESSING FROCTDURES 

C II. S. CEOLDGICAL SURVEY, ICNVER. COLORADO 

C EiRANCM OF FETROFHYSICS AND REMOTE SENSING 

C DON L. 8AUATZKY 

C 

c 

c 

C CE0MX4 GENERATES A RECTIFIED IMAGE FILE FROM A DISTORTED IMAGE 
r FILE AND FROM COEFFICIENTS FOR RECTIFICATION. R AND S, DETERMINED FDR AN 
C AFFINE TRANSFORMATION OF THE DISTORTED FILE. 

C INPUT FILE STRUCTURE CONSISTS OF A HEADER RECORD CONTAINING TWO 
C INTEGERS FOR LINE LENGTH. LENREC. IN PIXELS AND NUME'ER OF LINES. 

C NORECS. NORECS NUMDER OF DATA RECORDS FOLLOW, EACH RECORD CONTAINING LEfJREC 
C DYTES OF 8-DIT DATA. OUTPUT FILE HEADER RECORD CONTAINS TWO INTEGERS OF LINE 
C LF-NGTH. NPXDLIT, AND NUMBER OF RECORDS. LMAX. INPUT PARAMETERS. IPIXIN AND fJPI 
C XIN, ALLOW TAKING A SUP'En' OF THE HJPLtT FILE. OUTPUT PARAMETERS LMAX. MNP i X 
C , AND MXPIX ARE SELECTED ON THE LINE LENGTH OF THE INPUT FILE AND DEGREE 
C OF ROTATION REOUIRED FGR RECTIFICATION. SECTIONS OF THE OUTPUT FILE OF LENGTH 
C LMAX AND CONTAINING PIXELS MNFIX TO MXPIX ARE GENERATED PY ONE DR MORE ITER 
C ATIONS OF THIS PROGRAM. SECTIONS ARE CONCATENATED IN SUDSEOUENT PROCESSING. 

C SECTIONING THE OUTPUT FILE IS DONE IN RESPONSE TO THE MESSAGE*. "INBUF ARRAY 
C TOO SMALL." 

C 

C DECLARATIONS 

REAL R(3).SO) 

LOG I CAL* 1 1 NBUF ( 200000 ) . OUTGUF ( 3000 > 

INTEGER FCDINISSI.FCBOUTOS) 

C 

C SET PARAMETERS 

WRITE<6.99) 

99 FORMAT! IX, 'ENTER PIXEL/LINE COEFFICIENTS '«) 

READ(S.96) R.S 

WRITEtt.yS) 

98 FORMAT! IX. 'ENTER INPUT F RST PIXEL. NO. PIXELS') 

READ(5,9&) IPIXIN. NPIXU 
WRITE(G.97) 

97 FORMAT! IX. 'ENTER MAX. O IT PUT LINES, MIN/MAX PIXEL 'D 

READ!5,96) LMAX.MUPlX.riPXX 
96 FORMAT !G1 6.0) 

C 

C OPEN DATA FILE TO Tf^'iNSFORM 

READ!E) LENREC. NORECS 

C SETUP WORK ARRAY 

NRECS=MI N ( 200000/LENREC. NORECS) 

MXLINE=NRECS 

MNLINE=1 

J1-LENREC«M0D! l.NRECS)*! 

J2=J1-*LENREC-1 
DO 90 I=1,NRECS 

90 READ!S) (INBUF!J),J=J1 J2) 

C 

C.....OFEN OUTPUT DATA FILE 
NPX0UT=MXPIX-t1NPIX+l 

1F!NFXOUT.LE.SOOO. AND, NPIXIN.LE. 3000)00 TO IlO 

100 STOP ' 7DATA FILES EXCEED BUFFER WIDTH.' 

110 WRITE! 9) NPXOUT.LMAX 

C.-..., RE AD/ WRITE LOOP 

DO 210 LINE=1,LMAX 
DO 200 1PX=MNP1X. MXPIX 



INF»X-R<1M(IPX) 4 R(2>«<LINEI « R(3> 

IFMNPX.L7. 1.PR. Jr;PX.07.WPJXlN) 7HEN 
PUT PUr ( 1 P X - MNP I X ♦ I » • . F ALEE . 

ELSE 

INLlNE>S(n«()FX> 4 S<3)«(L1NE> ♦ St3> 

IFdNLlNE.LT. I.OR. INLINE. 07. NOREC3) THEN 
OUT PUr- < 1 PX -MNP 1 X 4 1 ) • . F AL&E . 

ELSE 

C CHECK LIFT FOR ECANLINE IN WORK ARRAY 

IF IINLINE.OE. online. AND. INL INE. LE. MXLINE ) 00TPJ20 

C ELSE HEAD FCANLINE INTO LIST AND WORK ARRAY 

1FUNLINE.lt. MNLJNE> STOPUNBUF ARRAY TOO SHALL!!' 

DO ilO l«nXLlNE4l . INLINE 
IREC«HOD< 1 ,NRECS>*LER’REC 

115 READ(S.rvEC*INLlNE> UNDUFI J> , J^IREC^) . 1REC-*NPX1N) 
rtXLlNE- INLINE 
NNL I NE » HXL 1 NE-NRECS4 1 

120 OUTD'JF ( 1PX-MNP1X4 I) lNt<UF < LENREC«r.OD ( INLl NE. NRECS ) + 1 NPX ) 
END IF 
ENDIF 

200 CONTINUE 

210 WRITE<9> <OUTBUF(J). J»1.NPX0UT) 

C 

C FINIS 

300 STOP 
ENO 


PFO'l-RAn HCMTir* 

C THIS FFOC*RAn CC'Mf 'JTE' A RELATIVE THERMAL INERTIA AND TEMR DIFF 
C IMAOE FILES FROM HCMM DATA. 

C 

DIMENS I ONI J Cr-). I2C:r.) . I S' ( SE ) , I 4 <?E. ) . IDIC'OOO) . ItlMl'CiOO). 
e- ii'?<iOi‘X'> 

DIMENtlON TEMP'C'E^^I.ALDd'Et'). 15<:<E) 

D 1 MEN : I ON I i- ( ' E > . 1 DUM I 200A ) , 1 DT < 2000 > 

DATA ii/:-*o. 1 . i.S' 0 «o/. i 2 /!*o, t. i.:' 0 »<i/. i:/:-*o. 1 . 1 , :-' 0 *o/ 

DATA 1 4/ It. 0.0. 1 . 1 .00*0/. 15/”?. 0. 0. 1 . 1 . 70*0/ . 1 < /0 . 0 . 0 . 1 . I . 20*0/ 
INTE'.'ER*!' IT <2000) 
c 

C OFEN NOON THERMAL FILE 
C 

WR1TE<« .RMT^M-" OPEN DAY THERMAL FILE")') 

CALL DI')’.IO (0. 12. IDl . 1 1 ) 

C 

C OPEN NIGHT THERMAL FILE 
r 

WR1TE<6.FMT»‘'<'" OPEN NIGHT THERMAL FILE")') 

CALL D1:)'I0 ( 0. 1 3 . I D2. 1 2 ) 

C 

C OPEN NOON ALDEDG FILE" 

0 

WR3TE<t.FMT*'(" OPEN NOON ALDEDO FILE")') 

CALL DITKIO < 0. 1 4 . 1 P? . 1 3 > 

C 

C OPEN thermal inertia ootput file 
r 

WRlTE(t.FMT»' ( " OPEN THERMAL INERTIA OUTPUT FILE")') 

Ul=! 1 (2) 

>;=I2(2> 

Y=1 :-<(2> 

W 1 = 1 1 ( 3 ) 

X1=I2(3) 

M = 1 2 < 3 ) 

MAXP=AM1N1 (N, X. Y) 

MAXL«AMIN1 (W1 , XI , Y1 ) 

I4(2)»=MAXP 

I4<3)*MAXL 

CALL DISf.IO <0.1?. 0.14) 

C 

C OPEN TEMP DIFF OUTPUT FILE 
C 

HRITE(6.FMT = ' < " OPEN TEMP DIFF OUTPUT FILE")') 

I?<2)»MAXP 

I5<3)=MAXL 

CALL DISKIO <0. 16.0, I?) 

C 

C OPEN DUMMY FILE 
C 

WRITE<6.FMT«'(" OPEN DUMMY FILE")') 

I6<2)<=MAXP 

I60)=MAXL 

CALL DIOKIO <0, 17,0. U ) 

C 

C INPUT AVERAGE PARMO 
C 

WRITE <6.FMT»'<" ENTER AVERAGE VALUER A.TD.TN")') 

ACCEPT * , AAVE. ATP. ATN 
TYPE • . AAVE. ATD, ATN 

WRITE <6.FMT«'<" enter LATITUDE AND GOLAR DEC IN DEC")') 




ACCEPT • .KL AM, DELTA 
TYPE • .KLAMr DELTA 
XLAM«rLAH*0.0»74f? 

DELTA-DELTA-0. 0! 74E.? 

£€•;*• 1 . / < C 0 ; ' • >.TM-DEL f A ) > 

0(m i ( j . - . 2*S OPT < L'EO > * 

C»-C»>*is^2.« 1P4*<C0?.< XLAf1)*0, •DELTA*? 1 N ( XLAM > ) 

&■ ( ATD- ATN ' • > D'.*0-C* AAVE 

TYPE • .00 

TYPE • .C 

TYPE ♦ .& 

P/MTE <6.FMT*M" P PEL • < P*C *ALD ) /DT " > ' > 

r. 

C COMPUTE TEMP AND ALD ^ALJPPATION FILES 


2 ALP<K)«<K-i; ■’0.00:'<«'2l57 

Cl*144il.V:;7 
. lT-91 

C3»-UP.2J:'-76 
DO i i LK»1 > O'E'i' 

1 1 temp<l>:)«ci'/<aloc.<ci/(l»;-i*-c:;)41 > ) 

c 

C MAJN P'EAD/WPJTE LOOP 
C 

DO 0 1«1.MAXL 

CALL D!-:nO <«', 12, IDl . H > 

CALL UNFAO'.l IDi , 11 (2) ) 

CALL DlSl'IO <P. I?, ID2, 12) 

C ALL UMP Af t ■ < 1 D2 , 1 2 ( 2 ) ) 

CALL DlSl lO <®. M, IDO, 13) 

CALL UN<='ACK< 1D3, i:-'.(2) 5 
C 

DO 4 I'.el.MAKP 
NUM1*1B1 (K)*l 
NUM2=ID20:)*1 
NL'M?=IP3(K)*l 
TD=TEMP(NUM1 ) 

TN«TEhP(NUM2) 

AA=ALD(NUM:-) 

IDUM<K)«2?5 

E-TD-TN 

X»:<AA-AAVE)*!00 

Y=ATD-TD 

IF <X.0E.2.AND.V.0E. iro DO TO 21 
IF lTN.LE.26tl> CO TO 21 
2 = -S 

IF (S.LE.O. AND. Z.LE. 10) CO TO 22 
IF (S.LE.O) 00 TO 21 

1 DT < K ) ■$ 

IF < lT(K).ECl.35.E-0) 17(K)*?5C0 
IF (IT(K).EQ. 1000) 1T<K)«?90 
00 TO 4 

21 IT<K)'1000 
IDT 0<)«i!?5 
1DUM(K)»0 

OO TO 4 

22 IT(K)«30?.0 
IDT<K)-0 
IDUM(K)e255 

4 CONTINUE 

CALL Dl'iKlO < 10. 1E>, IT, 14) 


CftLL PA*:» tlt'T,l5{2>} 

CALL ( io. u . iDTi ir.> 

CALL PACK( IDUM. 16<2) > 

CALL MfKlO no. 17. JI'UM, !<>) 
y CONTlNUt 

c 

C CLO-E FILES 

r 

CALI DJinO << . II'.C*. 1 1 > 

CALL C'lSKlO <S. 1 : 1 . 0 . 12) 

CALL MS) 1C <« . M.O. 12) 

CALL DISKIC <6. 1^.0. M) 

CALL MS) 1C << . n .0.1^) 

CALL MS)\C <0.17.0.14.) 

STOP 

END 
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JO t PKOFLt Ptofilin.^ r«>utinet. J/2/81 
20 1 Pcojtam l«n^u« 4 «r i>t tatvnoca Dabis. 

20 I Program wiitlcn by i(.i««t»on 

4C OPllON bASC ) 

SU OlH Tuay(lSOU) ,Tnilc(lS0U} « Clev (1 SuO ) ,x (1 &00) , y (1 SOO) ,Dum(l 500) »0* r v ( 1 ‘>00) .] nert(15C 
- OU) ,01(1 • Tcpo( ) 5u0) ,Mean(250u) 

»0 D 1 .I H(2uuu) ,Z(2000) ,X2 (2000),Topa(2U0U),S$te0),A5(5) t4lt 

70 IN'IOULH Dn(40uu) 

60 PRlN'rLR IS 26 

dU Llvv coir«0 I Inuiajizc c2ev«tion cortcction inaicator 

200 CnaoTfc'O I Iniv talize lot piot>2e plottin>j 

2 20 INPUT *k:ntet prctiJf: le. A , B” . Prot 1 2c$ 1 2 , 2 ) 

22o proi i2i^>Proii2eSlPiot i2c$b* '* 

2 jo iNPu'l ’l-nier area name ie . PRB* ,Atca$ 

2 40 Ki2e__Bp6«At ea$4Ptoi i2e$ ( 1 , 2 ) 

250 Function aeyr.: I 

2b0 On"XEY iO goto lopo_oen 

2 70 ON KCy 1 2 GOTO 2r.a9e 9en 

2e0 ON KEY |2 GOTO Proi plot 

2 5U ON KEY 1 3 GOTO FratUrc 

2o0 ON KEY 14 GOTO Fi2e_cdit 

22 0 ON KEY 1 5 GOTO CrosB_^cor 

22o ON KEY lb GOTO Cl Ota~plot_Sn 9 1 

230 ON key |7 goto CroBs“plot“i> ^2t 

240 Oi« nEY i6 goto Topo Tesample 

25u ON KEY 1 9 GOTO Elev”correct 

2bu UN KEY I 20 GOTO Repot t__p2ot 

27; 2F Eiev coir»I THEN Pr 01 1 le$ » Pr or i lejb * elev cotr* 

2b0 MASS STORAGE 2S "tC* 

290 EXIT GRAPHICS 
300 PRINT PAGE 

310 PRINT • KEY lO Generate topo binary tile from PERKIN output* 

320 PRINT * KEY |2 Generate irage binary tile* 

330 .PRINT * KEY |2 Plot a prorile* 

340 PRINT * KEY |3 Select a teature tor matching protiles* 

350 PRINT “ KEY 14 File edit* 

360 PRINT * KEY 1 5 Cross correlate two images* 

370 PRINT * KEY |6 Cross plot two data f i les (r epor t ) * 

360 PRINT * KEY 47 Cross plot multiple data tiles* 

390 PRINT * KEY 46 Resarr.ple the topo data to match images* 

400 PRINT * KEY |9 Determine elevation correction to temp data* 

410 PRINT * KEY 110 Profile plot t ing ( tepor t) * 

420 PRINT * MASS STORAGE SET TO ;C* 

430 DISP "Select option* 

440 GOTO 4 30 

450 1 

460 I 

470 I 

460 Topo gen: t To generate topo binary tile from digitized records. 

490 OVERLAP 

500 BUFFER ll 

510 PRINT PAGE 

520 INPUT *Enter number of input tiles from D1GIND*,K 

530 File5-*ToppBB* I * * 

540 D15P *Generating vopo file.' 

550 I-O 

560 FOR N»1 TO M 

570 ON END 11 GOTO 640 

560 ASFiCN 41 TO FileS 


n 


»tO READ IliK 

too RCOIM R(K) 

ilO READ lltCurvcS 

t20 READ lltXlot^cr .Xupper 

»30 READ ll:R(*) 

t40 ASSIGN II TO • 

«S0 FOR L-1 TO K 

itO F( 1*2) "RCL) * 13* .0234 I Convert iron ft to la, 

(70 NEXT L 

(00 X*1*K 

(90 NEXT N 

700 REWIND •;T13" 

710 REOIM Z(l*2) 

720 Z(l)-1 

730 Z(2)*.1387S I Xm . Aetuinex 1/40 * kpecing on 1:230,000 map 

740 Fileout5**Topc"lP lol ile$ |1 , 2) 

730 FCREATE Fileout$,52 

760 FPRINT Fileout3,2 (*) 

770 DJSP ■FinjBhed* 

780 MT 300 

790 SERIAL 

800 GOTO Function keys 
810 I 

820 I 

8 30 t 

840 Image gen: I To generate image binary files. 

850 OVEf^LAP 

860 SUPPER II 

870 PRINT PACE 

880 1 Regular Perkin image files 

890 I 

900 LINPUT 'Enter file type ie Tday ,Tn i te , Ref 1 * ,F i lex 3 
910 DISP 'Generating image tile.' 

920 XS»Prorile5(l ,l)k';Tl5' 

930 IF Filex5»'Tday' THEN F3»'TD'lX$ 

940 IF Filex3»'Tnite' THEN F3-'TN'iX5 

950 IF FilexJ-'Ref 1' THEN F6-'RD"4X5 

960 INPUT 'Enter starting pixel and line',Ps,Ls 

970 INPUT 'Enter ending pixel and line',Pn,Ln 

980 INPUT 'Enter number of oatr blocks ',No blocks 

990 Cl-(Ls-Ln)/(Ps-Pn) 

1000 J>1 I Profile element counter 

1010 1-1 t Block no 

1020 FOR 1-1 TO No Blocxs I INSERT FAST ALGORITHM HERE. 

10 30 CAT TO AS$*);fS4VAL5(l) 

1040 IF LEN(AS(l))O0 THEN GOTO 1090 
1050 REWIND ':T15' 

1060 PRINT 'F34VAL3(1) not on tape. Insert correct tape and COST' 

1070 P^USE 
1080 GOTO 1030 

1090 ASSIGN fl TO F54VALS(1) I 

1)00 READ «1;S$ 

1110 N-POS(S3, 'Lines') 

1120 Lines«VAL(S$(N+5,N47J) 
li30 M-POS(S5.'Pi*els') 

1140 Pixels»VAL(S$ |M46,M*B; ) 

1150 PRINT 'Lines, Pixels-';Lincs, Pixels 
1160 REDIm Dn(Pixels*Lines) 

1170 READ il;Dn(*) 

1180 FOR K-1 TO Pixels 




CRlGmAl PAGE IS 
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llvo vy»CJ*(K-)>4i 
12U0 U*lNT(yy) 

1210 ir yy-Ll>.& THEN L1«L1«) 

1220 (l.l'lJ*PlJiclk4K) 

4230 I PRINT J,Z(J) I Ukt {or editing. 

1240 J«J4l I increment profile counter 

12S0 NEXT K 

1240 ASSIGN el TO • 

1270 NUT 1 
13U0 REWIND *:T1S* 

1390 JbJ- 1 I Reset pxxel counter 
13u0 REOIM Z(J) 

1310 Kl*]4421.S07 I Temp celib Dn to Temp 

1320 K2»12S1.1S91 

1330 K3--116. 21376 

1340 Disf5UR( (Pn-Ps) *24(Ln-Ls)*2)*.4771 7/J I Xm spprox units 

13S0 REDIM Z(J) ,X(J) 

1360 IF (Filex5»*Tdsy*) OR (Filex$»*Tnite*) THEN GOTO Temp csl 

1370 IF Filex$«"Retl" THEN GOTO Pefl csl 

1360 Temp csl: FOR 1»1 TO J 

1390 Z(1F-K2/1.0C(K1/(Z(1 )-K3)4l ) 

1400 IF Z(I)<260 THEN Z(1)-26U 

1410 NEXT 1 

1420 GOTO Bin store 

1430 Be£l cs1:"F0R X-J TO J 

1 440 Z( l7-zm/266«100 I t refl 

1460 NEXT I 

1460 GOSUB Bin^stote 

1470 GOTO Function_^Xeys 

1460 Bin_store: I "" 

1490 DISP "Binsry store* 

1600 Dum.(l)>J i No^vslues 

1610 Dum(2)»Dist I Dx incremient in Xm. 

1520 FOR 1-3 TO J42 

1630 Oum(I)-Z(l-2) t Normsl ststemenr. 

1640 NEXT 1 

1550 REDIM Dum,(j42) 

1560 Filenames “Fl lex 5 t l,4)6ProfileSUt2] 

1570 FCREATE Filenames, 52 
1580 FPRINT Filenames, Dum(*> 

1590 PRINT PAGE 
1600 DISP ‘Finished* 

1610 SERIAL 
1620 RETURN 
1630 I 

1640 I 

1660 I 

1660 Prof plot: I On CRT or at 1:1,000,000 on 9872S. 

'1670 INPUT ‘Select plotter option; 0 CRT , 1 9872S* , Plotter 

1680 PRINT PAGE 

1690 COSUB Filu read 

1700 GOSUB File'plot 

1710 IF Plottcr-0 THEN EXIT GRAPHICS 

1720 If V-1 THEN GOTO 1690 

1730 GOTO Function_ke/s 

1740 File read: INPUT "Enter file type ie Tday ,Tn i ic , Ref 1 , Ele v , De rv , Dt , I ner t ,Mpan* ,F i le 5 
1760 REDTk XllSOO) ,y (1600) , Z (1600) 

1760 Filename S-File$[l , 4 ] iPr of i leS | 1 , 2] 

1770 IF (Fil€5-*Tday*) OR (F i leS -*Tn i te * ) THEN GOTO 2080 
1780 IF (FileS-'Refl*) OR (F i leS- *Ele v* ) THEN GOTO 20 80 


Cjo 



1790 IF rilc»«*DcrV TilCN C01O 20ti0 
laoo IF r»lcl>*Mean* THEN GOTO 1920 

UlO IF (rilcl-'Ot*) OR (nle|-*lntrf) THCI) GOTO Ot (oim 
lfe20 Dl 10(m: FRCAD *To*y *9Pc Of i lt$ 1 1 , 2) , X (• ) 

1930 J7x(l) 

1940 Dl»l>X(2) 

IfbO FXrAD ”Tnlt*4Frol il*(|l, 2| ,y (■ ) 

1>»6U R£DIM X(J«2) ,y (J42) |Z (J«2) 

1970 If (F>ie9»"Dt‘) OR (Fil#|«*lnetf) THEN MAT 2»X-y 
1100 IF r»l«$»"Mt*r* THEN MAT 2-X^y 
1090 IF F>lC$>*MC*n* THEN MAT 2>Z*uS) 

1900 IF (F»ic9»"Dl*) OR I F i leS • ' Me«n’ ) THEN GOTO 21)0 
1910 Inert (otmt FREAD ” R*1 l*iProtnc$ (1 , 2 ) , X ( * ) 

1920 IF A7e*i»*PRD* THEN GOTO 19b0 

1930 PRINT "Neeo to reset elpna.Deta corf i at inert (orn^* 
1940 PAUSE 

19S0 Alpna«9226!).0 I Powcer River Basin Parameters. 

19()u Beta*-61744 

1970 F.AT X»X*1.01) I Convert to fractions 

)9b0 FOR L-2 TD J-*2 

1990 IF Z(L)>0 THEN GOTO 2020 

2000 Z(t)«25b0 

2vl0 GOTO 2060 

2020 2(L) -(Alpha+eeta*X(L) )/Z{L) 

2030 IF Z(M<bO0 THEN 2(L)-4bO 

2040 IF Z(M>3000 THEN 2(L)-3000 

2060 DaTiiZ) >Dist I Da increment i, Kin, 

2060 next L 

2070 GOTO 2110 

20e0 FREAD Filenames ,2 (*) 

2090 J>2(1) 

21C0 Di£t-Z(2} 

2110 FOR I«1 TO J 
2120 Z(1)-Z(142) 

2130 X(l) »(I-1 )*DiSt 

2140 NEXT 1 

2)50 REDIM X ( J) ,2 fJ ) 

2160 RETURN 
2170 1 

2160 1 

2190 I 

2200 File plot; SERIAL I IF V-1 THEN REPEAT 

2210 IF Enable-1 THEN GOTO File cont 

2220 IF Plottei^O THEN PLOTTER Ts 1 3 , 'GRAPH 1 CS* 

2230 IF Plotter-0 THEN GRAPHICS 

2240 IF ?lotter-0 THEN GOTO File cont 

2250 IF V-0 THEN PLOTTER IS 7,5,"’9e72A" I First time 

2260 File cont: BEEP 

2270 Dz-‘5 

2280 IF (FileStl ,1 )-■£• ) OR (F i le$ |1 , 2 1 - 'To' ) THEN Da-100 
2290 IF FilcSll ,2)-"De" THEN :'?'10 
2300 IF Fj leSIl ,21-*ln* THEN Dz-500 
2310 yS-FileS 

2320 IF (ySll,ll-*T*) OR ( y $ 1 1 , 1 J -*D* ) THEN yS-ySkMK)’ 

2330 IF ySll,ll-'M* THEN yS-'Mean temp(K)“ 

2340 IF ySllflJ-’E' THEN y5-yS6*(n>>* 

2350 IF y$ll,ll-*R" THEN yS-ySfc'm* 

2360 IF yS|l,2J-’De’ THEN yS-"Gr ad ient (m/Km) • 

2370 IF y5(l,2j-’Df THEN y$»"Temp dlff(K)* 

2360 IF ysiljlJ-'l" THEN y S-' I ner t i a (TI U) * 




•i 



2310 xi«*Oifttanct (km) * 

2400 HAT SEARCH Z (• ) «HAX > ZiMN 
2410 MAT SEARCH Z (• ) ,H1 N > Zmi n 
2430 Zmin>INT(Zmin/Dt) 'Dt 
2430 Zmaa-(INT(Zmax/Di)«) )*Dt 
■ 2440 DsB*(Zm«ii-Zir in) * . 12 5/B . S 
2430 OKft-X(J)*.12S/S.S 
2460 rOR 1>1 TO J 
2470 X(l)-(l-l)*Di»t 

2480 NEXT 1 

3490 IF Plottrr«0 THEN GOTO 2670 
3300 DEC 

2310 LIMIT 0<400, 0,283 
2320 LORO 5 
2330 CSIZE 2.73 

2340 Xl«(l. 3*23. 4-12. 5)*100/283 

3330 X2-Xl*X(J)* 100/283 I Ptopet kcalt 

2360 Y1>(1. 3*23. 4^14) *100/283 

2370 y2>yl4)30*100/263 

2380 MOVE (X14X2) *.3, yi>. 3*23. 4*100/283 

2390 LABEL X$ 

2600 LOIR 90 

2610 MOVE X1-. 9*25. 4*100/285, (yi*y2)*. 5 
2620 LABEL V8 
2630 LDIR 0 

2640 LOCATE XI ,X2,Yl ,Y2 
2630 FRAME 
2660 LORC 5 

2670 SCALE 0 ,X(J) ,Zmn,Zmax 
2680 OUTPUT ?05;“VS2j" 

2690 MOVE 0,Zmax 

2700 DRAW X (J) .Zr.ax 

2710 dram X(0) ,Zmin 

2720 AXES lO.Dz.O ,Zmin,l0,2 

2730 FOR 1"1 TO J 

2740 IF 1-1 THEN MOVE X(l),2(l) 

2730 DRAN X (1) ,Z(1 ) 

2760 NEXT 1 

2770 OUTPUT 703; *VN;’ 

2760 CSIZE 2 
2790 LORC 6 

2800 FOR ysBZmin TO Zmax STEP Dz 
2610 HOVE -Dxs,ya 
2820 LABEL VAL$(ys) 

2830 NEXT YS 
2840 LORC 6 

2630 FOR N-0 TO X{J) STEP 20 
2860 MOVE W,Zirin-Dzs 
2870 LABEL VAL$(W) 

2660 NEXT W 
2690 LORC 3 

2900 MOVE X (J) * . 3,Zmax-Dzs*2 

2910 LABEL ’Profile * 6Prol i le $ (1 , 2 ) 6 * ' ”6* ”6Area$6* 
2920 IF Plotter>0 THEN LABEL "CONT to continue* 

2930 IF Plotter«0 THEN PAUSE 
2940 PEN 0 

2930 IF Plottet-0 THEN EXIT GRAPHICS 

3960 IF Plotirr-0 THEN GOTO 2990 

2970 IF X(J)*23.4>7 THEN OUTPUT 703;"AF;* 

2980 IF X(0)*25.4<-7 THEN OUTPUT 705;*AH;* 




Area* 



29i^l INPUT *Lnttf 0 to UNO *nti 1 to plot •notntc proiil**(V 

3000 ir </•) THIN nCTUKN 

30)0 ir Plottt(*0 THEN C010 End) 

3030 IF )MJ)«25.4>7 THEN OUTPUT 70i.j’Arj* 

3030 IF X<J)*35.4<»7 THEN OUTPUT 70bj*AMj* 

•3040 Enul : GOTO Function fttys 
30 bO I 

30t»0 I 

30 70 I 

30o0 Featuiei GCSUIb )il«^i»id 

3090 Ftat: I UPUT *Ent«{ Tc«turt *nU «ppr on pixel velue" tK$ , J P 

3)00 JJ-Jt-12 

3))0 IK J)<) THEN J)>) 

312u J2«Jp*)2 

3)30 IP J2>J THEN J2«J 

314 0 Sp-0 

3) bO INPUT *Jte CONT to list «na 1 to plot”tSp 

31o>J IF 5f«0 TnCN GOTO Llfct 

3170 W,DIK A2 (J2-J)-») ) ,£>un (J2-Jl4) ) 

3)6u t*0 

3)ti*U KOR K-Ul TO J2 
320u L-l/»l 
3210 X2(L)>H 

3220 Du.T.(L) ■1(K) 

3230 NEXT K 4 

3240 Lr«x>L 

32b0 MAT SEARCh Duxr.l* ) ,MAX| Rirex 
32b0 ^5AT SEARCH Durr.(» ) ,MINj Rmi n 
32 70 Dr« R.’TUx-R/Tin 

32e0 Feature plot; PLOTTER IS .’ ’.•GRAPHICS" 

3290 GRAPnli^S 

3300 LOCATE 1U,)00,20,90 

3330 SCALE J3 ,J2 .RiTin.Rra* 

3320 GRID 1 , (R.T,ax~R«.in) • . 1 ,J1 , Rmln 
3330 FOR L«1 TO Lir.ax 
3 340 IF L-3 THEN MOVE X2 (L) ,Du.Tt(L) 

3350 DKAn X2 (L) .Our.(L) 

33t0 NEXT L 
3370 LORG 6 

33oU FOR L-1 TO Lirax SlEP 2 

3390 MOVE X2(L> .Rffin- (Rrax-Rnin)«.0b 

3400 LABEL 0AL$(X2(L)) 

3410 NEXT L , 

3420 LORG 4 

3430 mCvE ( Jl+J2)».b,RiT,in-.04«Dr 
3440 LABEL Flle$ 

3450 MOVE Jp,RiT,ax'*’Dr*.02 
3460 LABEL Kb 
3470 LORG 1 

3490 MOVE X2 (1) . Rrr^xfDr* .02 
3490 LABEL *CONT to continue* 

3bo0 PAUSE 
3510 PEN 0 
3520 EXIT GRAPHICS 
3530 GOTO 3630 

3540 List: PRINT SPA(IO) , "Feature ";K$ 

3550 Min-)E99 

3560 FOR 1-Jl TO J2 

3570 Min*MlN(Z(I) .Min) 

3580 NEXT 1 




Or QUAUTV 



3»t0 rOR 1>J) TO J2 

Skoo ir ^moMin TitcN MiijT i.z(n 

3ti)0 ir k{n»Mn THEN RRI.^T * 

3t3i/ KtXT 1 

)»J0 iNfUT *Entir 0 to titminc otnti iratutct ; 1 to tern ir>*ic* iT 

3<>40 ir T*0 THLt. GOTO rixtutt 

3»S0 ir T*i THLh GOTO tunctton xryv 

SkiO I 

3*70 I ......................... 

3*»0 I 

3*^0 C( 0 »» cort I FROrt3 To ciott coifvlate pioltlinj data tor U'*)**. 

3700 I nclnodi Snitt mean* to 0 and variance* to 1 
3710 I roim cro** cotiel *un tor vafym 9 uelay*. 

3720 I 7a>^e oata in lOO pixel *et*( See Airax) 

3730 I-RIUT i-tsQl 

3740 INPJT •tr.iet iunctionl I* fjay , Tnl te , Red 

37b0 lUFJ'i "Enter lunction3 ie Toay , Tni te . Ref 1 ...."iVS 

37c0 OISP "Ctot* cor ellat ion .* 

3770 FRLAD V S 1 1 , 4 ) *P lol i le$ 1 1 , 2 ) , R (* ) 

37*0 IP (a2(1)«R(1)1 MD (X2(2)-R(2)1 THEN GOTO 3B30 
37S0 PRltir "Pile* are not aame type - recnecx* 

3eo0 PRINT " Een, tn»" j X { 1 ) , X (2 ) , "D i *t«" ; R (1 ) , R( 2) 

3elu INPUi' "Do you wi*n to continue anyway 7",Anz3 
3e2U IF Anzillil 1»*N" THEN STOP 
Je30 J>x20) 

3e40 r (2) 

3dbC I ax; Kr^ax^u 010 100 
3ecU J ‘lou 
3o7o PnINTEK IS 0 
3etu Pts liil 

3oi»u FMn‘ 1 “ Eist OX correlation for "sX5;" ano "jyS}" P rol i ie" ; Proi il# 5 

5iv0 PR INI 

3il0 FOR h»l TO nrax 

3»<o 11 • (H-1 ) • lo04l 

3S3J 12-11+J-l 

3»40 PRINT "Frc/iT pixel " t v/ ALS ( Il ) * " to "*0X^5(121 
39E0 fOr 1-11 TO 12 
3960 X(l-ll+l) ■a2(W2> 

3970 V(l-n + l)-R(M2) 

39oO NEXT 1 

3990 REDIK X(J) ,y (J ) ,Dua(J) ,2(J) 

4uu0 K,canx-SUM(X)/J I Transiorffi to Mean-0 Dev-1 

4010 Meany-SUIMY) /J 

4020 MAT X-(-Htanx)+X 

403o HAT y-(-Meany)+y 

4U40 MAP Dur.-X.X 

40b0 Vat-SU' 1 (Dutr)/J 

40o0 HAT X-(l/S-wR(Var) )"X 

4070 MAT Durj«y.y 

40e0 Var-bUM(Dutr.)/J 

4090 HAT y- (l/SwR(Var) )-y 

4100 PRINT 

4110 PRINT " Delay *;X5!".";y5 

4120 FOR Delay--S TO 5 

4130 Suira-0 

4140 FOR 1-fc TO J-5 

4190 D»I4Delay 

41o0 SURi-Afc5(X(l)*y(D) ItSuml 

4170 Next: NEXT 1 

4180 PRINT " " {Delay ,Suml 


9 ^ 



«1I^U 
420U 
4210 
4220 
423u 
*4240 
42 SO 
4240 
427U 
42»0 
42^0 
4300 
4310 
4320 
4 3 30 
4340 
43S0 
4340 
4370 
4340 
43»0 
4400 
4410 
4420 
4430 
4440 
44S0 
4 460 
4470 
4440 
4490 
4S00 
4S10 
4S2G 
453> 
4540 
4550 
45oC 
4370 
45d0 
4590 
4tu0 
4610 
4b20 
4 6 30 
4640 
4b50 
4660 
4670 
4640 
4690 
47u0 
4710 
47?0 
4730 
4740 
4750 
4760 
4 7 70 
4 7o0 


NLXT Dtlay 
NLXT K 

rKlNT£4 IS U 
C070 runcuon 4*yt 
I “ 

I 

I 


Cl okK_plot_tn| It I 
Plot“counT« 1 

01SP~‘''CiOES plotting > pinjlt.* 

lUPin *StUci plotur opttoni 0 CRT • 1 9472S* .P lottcc 
IF Plcttt{«J TMCN PLOTTER IS 13,'CMPHICS* 

IF Plotltf-l THEN PLOTTER IS 7,S.*9S72A* 

IF Plcttt(>l THEN LIMIT 0.400,0,265 
CfOkE conn INPUT ‘Enter filcnamt lor X«k>i le Tday ,Tnttf . . . 
input ‘Enter rtlcn«n<* tor yaiiik x» Toay.Tnitc. . ,y S 


Fk»-xS|1,2| 

Fy4-ySll,2| 

IF FkS|1,2)-‘To‘ then 
IF FySl 1 ,2J»‘To‘ THEN 
IF F*$|l ,2)»‘Tn‘ THEN 
IF Fyj( l,2)«‘Tn‘ THEN 
IF Fk6|l,2l«‘R«‘ THEN 
IF FySl 1 ,2J»‘Re‘ THEN 
IF F»$l 1 ,2J»‘E1‘ THEN 
IF FySll,2J*‘El‘ THEN 
IF FK611,2J-‘De‘ THEN 
IF Fy6|l,2)«‘Oe‘ THEN 
IF (F*S|l,2)-‘ln‘) OR 
IF (Fkkll,2i>*Dt*) OR 
IF (F«J|1,2)»‘.M«‘) OR 
GOTO Plotting ketup 


FREAD ‘Td«y‘tPronicS|l ,21 ,X( *) 

FREAD *Tday‘6Ptoi ileSi 1 ,2] ,y (« 1 
FREAD ‘Tntt*6ProtUtS|l ,2| ,X(‘l 
FREAD ‘Tnit‘4Pronit$|l ,2) ,y (• 1 
FREAD ‘Rcn‘kProf ilcSll ,2) ,X(‘) 

FREAD ‘Rcll‘kPronic$(l ,21 ,y(* 1 
FREAD ‘Elcv‘4Prolilc6|l ,2) ,X(‘) 

FREAD ‘Elcv‘4Profilc$|l,21 ,y(‘) 

FREAD ‘Dcrv‘4Prolilc6|l,2) ,X(‘l 
FREAD ‘D»rv‘4Proril*S(l ,21 ,y (•) 
(FyHl,21-‘In‘l THEN GOTO lntrt_^read 
iFy$ |1 ,41b*Di‘) then GOTO Dt read 
(FyS 11 ,4l»‘Me‘) THEN GOTO Dt"read 


Inert rcao: I 


Alpr<«>5 2265.6 I Powder River Bakin Parairieterk. 
Bcta>>61744 


FREAD ‘Toay"tProrlle511,2] ,X2(‘l 
FREAD •Tnit‘kProtne6|l,21 ,£(•) 

FREAD ‘Retl*4Protilc5|l,21 ,Duni(*) 

REDIM Duir.(ROn (Z) } 

HAT Duir^Dun.* ( . 01) I Convert Retl to a fraction 


KAl X2-X2-2 
MAT X2*x2*(Beta) 

MAT X2>X2^ (Alpna) 

MAT X2«X2/Diuti 
FOR 1-1 TO ROw(X2) 

IF X2(l)<1000 THEN X2(l)-1000 
IF X2(l)>3000 THEN X2(1)>30U0 
IF Fx6«‘ln* THEN MAT X-X2 
IF Fy5«‘ln‘ THEN MAT Y»X2 
COIC 4400 

Dt react FREAD ‘7cay ‘ 4Pr or i le$ (1 , 2 1 . X2 (‘ ) 

FREAD *7ni t'tProf ile$D , 2) ,Z '• ) 

IF (Fx5"’Df) OR (Fy^-’Pt") THEN HAT X2»X?'Z 
IF (Fx6«’Re‘) OR (FyS*''Me‘) THEN MAT X2»X?->2 

IF {Fx4»‘Jle‘) OR (FyS»‘Me“) THEN MAT X7 > 

IF (rx5»‘DV*) 0,' (Fx5-‘Me‘) THEN HAT X- 

IF (ry$-‘Ot‘) OB (Fy4-‘He‘) THEN MAT y-. 

GOTO 4500 

Plotting ketup: IF (Fx50‘ln‘) OR THEN J-X(ll 

IF (Fx5">‘ln‘l OR (Fx50‘Dt‘) Tn. X(2) 



ORKMML PAM M 
OP POORQUAUrt 


47VU 

4aU0 

4oJU 

4420 

4430 

4440 

44&U 

44t)0 

4r70 

4«60 

4b90 

4VuO 
4970 
4 920 
4 930 
4940 
4 99u 
49bU 
4970 
4 9e0 
4990 
&O&0 
30i0 
9u20 
50 3J 
50« 0 

50 50 
5060 
5070 
5040 
5090 
5100 
5110 
5120 

51 30 
514 0 
5150 
5160 
517y 
51b0 
5190 
5200 
5210 
5220 
523u 
5240 
5250 
526u 
5270 
52b0 
529C 
5300 
5310 
5320 

• 5330 

5340 
5350 
5360 
5370 
5380 


IF 

(Fa a 

• "1 r." 

) UR (Faa*"Dt") 

THEN J> 

•X(l) 

IF 

FxSa 

"He" 

TnCN J>Z(1) 



IF 

(Fa4 

• "I n" 

) OR (Fxa»"Dl") 

THEN Dl 

>St> Z( 2 ) 

IF 

Fa$> 

"Me" 

THEN Dist>Z(2) 



FOR 1>1 

TO 0 



X(1 

)-X( 

142) 




»(1 

)"V( 

1 42) 




NEXT 1 





REDIm X 

(J) 

(J) ,X2(J) ,Z (J) 



Dx* 

^lu 





Dy ' 

•10 





If 

Xbil 

.11" 

L" TnEN DxalOU 



IF 

yxii 

• 11* 

L" THEN DyolOO 



IF 

xbli 

JJ" 

1" THEN DX-5UU 



IF 

y6 11 

• U" 

1" THEN Dy»50U 



IF 

(X51 

1 . 2 ] 

"Tn" ) OR (X5 (1 , 

2J«-Dt") 

THEN Dx>5 

IF 

X$|l 

• 2J. 

Me" then Ox -5 



IF 

(V$| 

1«2) 

"rn") OR (y$n, 

2J»"Dt") 

THEN Dy-5 

IF 

y>ii 

,2J ■ 

Me" THEN Dy«5 



IF 

xali 

• IJ" 

R" TnEN Dx»5 



IF 

y$ii 

• IJ" 

R" THEN Dy«5 



IF 

Plot 

ter >0 

THEN LOCATE 2U 

,100,10, 

90 


COSOo foint (<lot 
csm 2 
LOKG 5 

MOvi: ( x>TtAx«Mun) * . 5 <y n.«x 
LAfiCL ”Protil* "iPtotile^ 

IF Plottcr>0 THEN GOTO 5210 

INPuT "Enter 0 >.o END and 1 to continue plots” »Sz 

IP Sz-0 THEN GOTO 5160 

IF Plot count«4 THEN GOTO Reset plot 

Plot count<"Plot count^l | increment plot counter 

COT &”C r o s s_^ CO n t " 

Reset plotr^OOTPUT 705;”AH;”J Next paje 
" OUTPUT 705;”AH;” 

Cross__plot >1 
GOTO Cross cont 
OUTPUT 705; "AHj” 

OUTPUT 705;”AH;” 

DISP "finished.” 

PEN 0 

GOTO Function xeys 
LORG 2 

MOVE Xmin ,ymax4yr 

LABEL "CONT to cont i nue ; DUMP GRAPHICS to copy" 

PAUSE 
PEN 0 

EXIT GRAPHICS 

INPUT "Enter 0 to continue and 1 to DIGITIZE plotted points"«Z8 

IF ZB«0 THEN GOTO Function xeys 

GRAPHICS 

POINTER (Xiran+xmax ) • .5, (ymin+ymax) *.5,2 
DIGITIZE Xc,yd 1 Position cursor 
MAT Z»X-Xd 

mat Z-Z.Z 
MAT x2-y-yo 
MAT X2-X2.X2 
PIAT 2-Z4X2 

MAT SEARCH Z (• ) ,MIN ; Ziri n 

HAT SEARCH Z (• ) ,LOC {Z (K) -Zmin) ; X 







S3lit> ^^VL (A.Rin«j(n«ii)*.&.yR4«4)rx 

S40& LAbLl "riMtl t*>vAL<(K)|* .U%t CCtiT to continut" 

»4Ji) fAUSC 
S42U Gcro &2bU 
&43J i 

5i4u 1 

&4M 1 

S 4 »u Cross plot inultt I Cross plots or Td«y tTni t» ,Topo etc. 2/20/bl 

S«7u J'BJnT fagI 

54eu "Svlsct plotter option: 0 CHT < 1 96723” , Plotter 

i 4 f>u D15P "Cross plotting > inultiplc.” 
iboO CU3U6 MUit_t ilc^reeu 
5510 GOTO Cont“inult” 
l-.ult t'lle Fesut i 
5530 Pb>7ror iTejil 1 ,2) 

5540 If Aren3«"PRC" THfN GOTO 5570 

555i; FRliiT *Neeo to revise P n^ipping perms et Pcoetl” 

53d0 PAL5C 

5570 Pcoetxt AlphsB52265.6 I PoBOcr River Basin Parameters. 

5so0 BetaB-bl744 

5690 FREAD "loay” 4P5 ,Toay I'*' ) 

5bu0 FREAO ” El e v" fcP5 ,TopO ( * ) 

5bl0 FKEAD **Tni t"iP$,Tnite I •) 

5620 FREAO "Oer v” SP5 , Derv { • ) 

5630 FFEAO "Ret I’fcPS ,Ren(«> 

5640 Oist- ro!^o(2) 

5b50 J«TopO(l) 

5660 FOR L>1 TO J 
5b70 X (L) •(L-1) ”DlSt 
5bd0 I'cay (E) BToay (L42) 

5b9o Tnite(L) -Tniie(L42) 

5700 TopO( L) BTopo{L42) 

5710 Derv(L) -Eer v{L42) 

5720 Dt(L) «roay (L)-Tnite(L) 

5 730 Mean(L) -(Taay <L) 4Tnite(L) )/2 

5740 fieil(L) -Rerl(E42) *.01 I Convert to tractions 

5750 Inert(L)-(Alpna+Beta*Re£l (L) )/Dt (L) 

5760 IF 1 r.ert (L) >2500 THEN 1 ner t (L) »2 950 
5770 If Inert(L) <1000 THEN Inert(L)-550 
5760 IF Retl(L)>.2 THEN Refl-,245 
5790 IF Oerv(L)>30 THEN Derv(L)*35 
5800 IF Detv(E)<-30 THEN Dctv(L)--35 
5610 NEXT L 

56 20 REDIM Tday(J) ,Tnite{J) ,Topo ( 0 ) ,X (0 ) , y ( 0 ) , Derv ( 0 ) ,D t ( 0) , T ner t ( J) , Ref 1 { 0) , Kean( 0 ) 
5830 RETURN 
5640 Cont mult: 1 

5650 IF Flotter-0 THEN PLOTTER IS 1 3, "GRAPHICS” 

5860 IF Plotter-1 TilEN PLOTTER IS 7,5,"9872A" 

5670 GRAPHICS 
5680 FOR K-1 TO 20 

5890 IF K KOD 4-1 THEN LOCATE 10,50,60,90 

5900 IF K MOO 4-2 THEN LOCATE 10,50,10,40 

5910 IF K MOD 4-3 THEN LOCATE 70,110,60,90 

5920 IF K MOD 4-0 THEN LOCATE 70,110,'10,40 

5930 ON K GOSUB PI ,P2, P3 , P4 , P5 ,P6 , P7 , P8 , P9 , PI 0 ,P1 1 , PI 2 , PI 3 , P 1 4 , P15 , P16 , Pi 7 , P 1 8 , PI 9 , P2 0 

5940 IF (K MOD 4-0) AND (Plottet-0) THEN PLOTTER IS 1 3 , “GRAPH ICS" 

5950 IF (K MOD 4-0) AND (Plotter-1) THEN OUTPUT 705j"EC;" I Enable cutter 

59b0 IF (K MOD 4-0) AND (Plotter-1) THEN 01- 'POT 705j"AF;" I Adv and cut 

5970 NEXT K 

5960 IF Plottet-1 THEN OUTPUT 705;"AHi" 1 Adv and cut 


»I90 

PLN 0 



itfoo 

EXIT CRAPHICS 



«U10 

D166 *riNlSHCO.” 



*02U 

GOTO Punctton kcyl 



‘ ftOlO 

PI I MAT X-Topo” 1 

Td«y vs Clrv 


«U4U 

MAT ¥«T(i*y 



*uso 

X6«"TopO* 




lf9*"T«i«y' 



»U70 

Dr *100 



ftUHU 

Dy«5 



6U90 

GOTO Cent 



fclOO 

P2; NAT X*Dcrv 

1 T<l«y vs Detv 


«no 

MAT ¥»TO*y 



ftl2U 

X6«’Derv* 



6Ua 

y$»“T0Ry" 



6140 

Dr«10 



61 bO 

Dy-5 



6lbU 

GOTO Cont 



6170 

P3: MAT X>X*tl 1 

Tdsy vs Retl 


6160 

N>T y»ro«y 



6190 

XJ-'Hefl* 



6200 

y6«*T0Ry* 



6210 

D»« .05 



6220 

Dy5 



6230 

GOTO Cont 



6240 

P6 : NAT X-Derw 1 

1 Tnttr vs Oerv 


6260 

NAT y»TniU* 



6260 

X$«“Dcr v" 



6270 

y5»"Tni le* 



6260 

Dx-10 



6290 

Dy«5 



6300 

GOTO coni 



6310 

P5: I<IAT X»ropc 

1 Tnite VE Topo 


6320 

NAT y»mite 



6330 

Xb«*TOpo* 



6340 

yS»*Tnj te* 



63b0 

Dx-lOO 



63t0 

Dy5 



6370 

GOTO Cont 



0 360 

P7: MT X»Aetl 1 

Tnile vs Rell 


6 390 

NAT y»rnite 



6400 

X6»’Re£l* 



6410 

V5»“Tni te“ 



6420 

Ox* .05 



64 30 

Dy«5 



6440 

GOTO Cent 



64 50 

P9 ; NAT X*Topo 

1 Ref 1 vs Tope 


6460 

MAT y-Re£l 



6470 

X,5»"Topo* 



6460 

y5*“Rell* 



64 90 

Dr-100 



6500 

Dy- .05 



6S10 

GOTO Cont 



6520 

Pll: NAT ji-lnert 

1 Refl vs Inert 


6530 

MAT Y-Reii 



6S4U 

X5-"Xnert" 



6550 

y 6-"Rctl* 



6560 

D«-500 



6570 

Dy-.05 



6560 

GOTO Com 

I? 

ORiCf'^IAL PAGE 15 
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Plot MAT A-oerv 

S«U0 

MAi y*6etl 

«»)0 

*6»"Derv* 

s*:r0 

y$«*Rell" 


0««10 

ftt40 

Dy«.05 

• *tsv 

GOTO Cent 

ftftbO 

PlSt MAT X»Ot 1 

6670 

MAT y •Inert 

6600 

Xj.-Dt" 

6600 

yC«»" 1 ner t* 

6700 

D«-5 

6710 

0y*500 

6720 

GOTO cent 

6 730 

Pi 3; hJiT X*To^l 

6740 

7JAT y»Inert 

6760 

X8-"Tof-o" 

6760 

y5»" Inerti*" 

6770 

Dm-100 

6 760 

Dy-500 

6790 

GOTO cent 

6600 

Fl4 1 MAT X-Oervi 

6 elO 

MAT y-lr.ett 

6620 

A6“"Der v" 

6830 

y6«*Iner tis" 

6840 

Ox-10 

6690 

Dy-500 

ob60 

GOTO Cent 

6870 

Pi 7.5 MAT X"Topol 

6 660 

MAT y-Dt 

6890 

X5«"Elev* 

6900 

y§.'*Dt" 

6910 

Dx-500 

6920 

Dy»5 

6930 

GOTO Cent 

6940 

Pie: MAT X«Dvrvl 

6950 

MAT y-Dt 

6960 

X5»"Derv" 

6970 

y5-“Dt“ 

6960 

Cx-10 

6990 

Dy«5 

7000 

GOTO cent 

7010 

PI 9: MAT X«Tnitel 

7020 

MAT y»Tasy 

7030 

X5«"Tnite" 

7040 

y$«"T6sy* 

7050 

Dx«5 

7060 

DylO 

7070 

GOTO Cent 

7060 

P4: 1 LABEL 

7090 

F8: i 

7100 

Pi 2: 1 

7110 

P16: 1 

7120 

P20: GOTO Plot^li 

7130 

Cents 1 “ 

7140 

point plot: 1 

7150 

MAT Search x(*) 

7160 

MAT SEARCH Y(«) 

7170 

MAT SEARCH X(*) 

7160 

MAT SEARCH y(*) 


I RtZl v» Dctv 


I nirt vs Dt 


Inert Vi Detv 


vs Elev 


Ot vfi Ders' 


7)*U i(i«tn>lNT(Kir.in/t>K) *Ds 
7300 X»i«K*(2KT(Xri«ii/DH)4n*Oii 
73)0 ytr.tn* Xiil (Yintn/Ly) *Oy 
7230 yir.«a'(ltiT(yin«>/Oy) 4) )*Oy 
03)0 Xt* (Xm*R-Miin) • .0& 

7340 yr» {yff,»R-ym»n> *.05 

72) 0 ir Plotltr>U THLN OOTO 7&40 
7240 LOKC b 

737u C;UC 2 

724(! if Plct count-) THCl. GOTO PU ) 

72 ; )F Pl&t“count«2 TrtCN GOTO Pll”2 
7300 )» P)ot"count-3 THLN GOTO Pit") 

73) 0 PU 4i Xl«f« . 5*2b.4 
7320 X)-Xfcett60 

7330 X2'X)43.7!>*25.4 

7340 y)>)!>0.7& 

7350 V2-254 

7300 GOTO 7520 
7370 PU )t X)-69.45 
?3tt0 X2>)65.) 

7350 y]-3d.l 

7400 Y2-)33.35 

74) 0 GOTO 7520 

7420 PU 2: Xl>69.ti5 
7433 X2«le5.) 

7440 y)-)5o.75 

7450 Y2-254 

74t)0 GOTO 7520 

7470 PU__3: Xie t -6 . 5* 2 5 . 4 

74d0 X) "Xset+eO . 65 

74»0 x2«X6et4) 65 . J 

7500 y)-3b.) 

7510 Y2-133.35 

7520 2z»100/(n*25.4) 

7530 LOCATE X) • 3Z , X 2* Zz , Y 1 ’ , Y 2* 42 

75<<0 SCALE Xm:n,Xr.«x ,Yr.in«Yn.ax 
7550 AXES Dx,Dy,Xmin (Ymin 
7560 LORC 6 

7570 NOOL (Xr, in4Xr.ax )/2# Yrr.in-Yr* 3, . 5 
75b0 FOR R-1 TO 2 
7590 IF R-1 T“iEIl 25-X5 
7600 IF R-2 V/IEN 25-Y$ 

76)0 IF 2S(1,1)-"M’ THEN 25»"Ke«n tei:.p(K)’ 
7620 IF 2S(1,1)»*T’ THEN 2$-256"{K)" 

7630 IF 25fl,2)-"Dt" THEN 2$»"Tcirp dUUK)* 
7640 IF 2511, 2) -"De" THEN 2 $• *Gr ad ien t (m/Xm) * 
7650 IF 2511,2)-*E1' THEN 25 • "Elev at ion (m) ’ 
7660 IF 25U»2)-*ln" THEN 25«" Inertia (TIU) * 
7b70 IF Z5[l,21»’Re" THEN 25-"Re£l(%)’ 

7680 IF B-) THEN X5’^i.5 

7690 IF R-2 THEN Y5-ZS 

7700 NEXT R 

7710 LABEL Xv 

7720 DEG 

7730 LOIR 90 

7740 LORC 4 

7750 MOVE Xmin-Xr*3. 5, (Ymin4Ytr,ax)/2 

7760 LABEL Y5 

7770 IDIR 0 

7780 f OR 1-1 TO J 


I(P0 



OMOMM. PMC M 
OF POOR QUALITY 


7710 MOVE K(I) .y (I ) 

7«00 LABEL 
7110 NlaT I 
7B20 CSUC 2.2 
7130 LORC fi 

7440 rOR A'Xmtn TO Xm«x STEP Ok 
7150 MOVE X.YiTUn-yr* .5 
7k60 LAECL 0AL$(X) 

7b70 NEXT X 
7660 LORC 8 

7b90 FOR YfVmin TO ViroK STEP Dy 
7900 R'.OVE xmin>xr« .5.VO 
7910 LA6LL VAL$(y«) 

7920 NEXT y« 

7930 CSUE 3.3 
7940 RETURN 
7950 Plci lAoel: I 
7900 SCALE 0 J ,0.1 
7970 LORC 1 
79B0 KOVE .2,. 75 

7990 IF AieAS>”PRB” THEN LABEL *Powtiei Rivtt Basin AtCA* 

6000 IF A(ea$>‘‘CP' THEN LABEL "Cabtzs Ptitta Aiea* 

8010 HOVE .5, .5 

8020 If Area5»"PRB* THEN LABEL 'Auj 20,1976’ 

6030 IF Axe»6»’CP" THEN LABEL ’Apt 3-4,1979’ 

e040 HOVE .5, .25 

6050 LABEL ’Ptof i It’SProl i le$ 

8060 IF Plottci>0 THEN DUMP GRAPHICS 
8070 RETURN 
6080 I 

6090 1 

6100 1 

6110 Top o_re sample: 1 To spline-smooth topo data and restore as Elev binary. 
6120 ! and Derv binary. 

6130 I Special correction to List and for topo ir.terp. 

8140 PRINT PAGE 

6150 DISP "Topo resample.’ 

8160 FPEAD "Toay’fcPromeSll,21 ,Tday<M 
6170 FKEAD ’Topo" *Proi i leS ( 1 , 2 1 , Tops {* > 

61 60 J>Tcay (1 ) 

6190 Dist»Tcay{2) 

8200 K»Tcp*(l) 

8210 Distz>Topx (2) 

8220 FOR I-l TO J 

8230 X (1) ■( 1-1 j’Dist I Image file 

8240 NEXT 1 

8250 FOR !■! TO X 

8260 Top*{ 1) -Topx(l42) 

8270 X2(I)><I-l>*Dl6tZ 

82d0 NEXT 1 

6290 t TO interpolate topo data 

8 300 DIM Xb{4 , 3) ,Xbtrn{3, 4) ,Durrj((3, 3) , Inv (3 , 3) , F ( 4 , 1 ) , Rum(3 , 1 ) , G(3 , 1 ’/ 

6310 MAT Xb«{l) 

6320 Xb(l, l)-( (0-XbO ,2) )-Xb(I,3) ) 

8330 Xb(3,l)«4 

6340 Xb(3,2)-2 

6350 Xb(4,l)-9 

6360 Xb(4,2)«3 

8370 MAT Xbtrn»TRH(Xb) 

6380 MAT Dumx«Xbtr n’Xb 


*310 lnv>lNV(Uuiw) 

11400 Dtl>DlSt/(2*Di»tS) I >0/2 

0410 L«1 

•430 lm«H>lNT(Difttl*K/Diftt) 

• 430 IF J>1»«K THEN I Truncatt all tilca to Jirax 

•440 rOF 1-2 TO J-1 

• 4b0 IF (X<l)>>X3(tM and (X(1 )<>X2 (L«1 ) ) THEN GOTO IiUctp 

• 440 L-L41 

•470 GOTO 1400 
4440 Inttip: F(1 >1) •Tops (L>1) 

• 400 r(2,n>TopR(L) 

•300 r(3,l)-Topx(L4l) 

•310 F(4,l)>TopK(L43) 

•320 NAT Xur>XDtrn«F 

• 330 NAT G*lnv*Run 

• 340 Zt>- (1-1 )*2*';tl>(L-2) 

• 330 Dur.)(l)>G(l .))/3* (Dcl‘243*Zb*2)4G(2, 1 )• Zb4G (3 ,1) 
a36U Dtrv(l)*2*G(l.l)*Zb4G(2,l) 

• 370 NexOt NEXT 1 

a3ei0 Dum(l )«Topx(l) I End points 

• 390 Oum(J) •Topx(J) 

•eOu Dcrv(l )-(Topx (2)-Topx (1) )/oiatx 
•olO Detvjj) >(Topx ( J) -Topx (J-1 ) )/Diatz 

• 620 XLGIM Dun.(J) ,Dcrv(J) ,Z(J) 

• 630 FAT Z*Cutn 

•64U n£OI.-) Topo( J42) iDerv (0v2) i Store output as Dinary tiles. 

••30 FOX 1*0 TO 1 STEP -1 



• ••0 Topo( 142) -Oumd ) 

6u7u Dervj 142) ■Deiv(l ) 

• 660 liEXT 1 
u69U Tcpo(l)*J 

• 700 Oerv(l) »0 
•710 Topo(2)-0ist 
11720 oerv(2)-Dist 

6730 KlJDlH TOfO(042) fDetv (J42) 

8740 DISP *5totino oinary tiles.* 

•730 FCKEATE *Elev" aPtOl ilej (1 , 2) , 52 
d7oU FCFXATE *Detv* a P rot i le $ ( 1 . 2 ) , 32 
6770 FPRIIIT “Elev'lProf ileS |1 , 2) ,Topo(* ) 

6760 FPRINT ”Dcr v* fcPtot i le$ |1 , 2 ) , De c v (• ) 

•790 DISP “FINISHED" 

6600 GOTO Function keys 
8610 1 
•820 I 
8630 1 

8640 File edit:! To edit binary tiles. 

6630 PRINT PAGE 

• 660 INPUT "Enter File to be wdited ie, Tday ,Tni te ,Ref 1« Topo. . , Fi lein$ 

• 67U Filenaire3«Filein$[l ,4]aProiilc$|l,2) 
odoO REDIM Dur.(130U) 

• 690 TREAD Fi lenan<e$ ,D uni( • ) 

•900 J*Dum(l) 

6910 Dist>Dum(2) 

•920 PRINT Filenames," Di5t»*;Dist 

• 930 PRINT LlMl) 

•940 PRINT •Dum(3)-*jDum(3),*Dum{"4VALS(J42)k")»";Dum(J42) 

•930 Edit teatute: PRINT "PAUSE AT Edit feature* 

•960 pau5e 

•970 Dum(l)-494 I 44***4«**4 SPECIAL MOO FOR PRB BB •••••*. 

•••0 0um(2)-.306 


\0t 



•f»0 FOR 1>3 TO 4»tl sp clud)* lot ' 
tdOU Duir(I)BDum(l42) 
fOlO NEXT 1 

• 020 RCDIH Duit.(496) 

•030 RLCP 

r040 PAUSE 

• 030 INPUT ”Uo ycu witr. to ft ttort tnit lilt 7*,A$ 

•U40 ir As|i#)|<>*Y' then cord mn (ttturt 

•070 PURSE t'ilcncDeS 
9000 FCRtATE F i lcn«ni>> , 32 
•uvO FPRIM r f Icntnu'S I Uum( •) 

•)0U GOTO function neyt 

•no j 

• 120 I 

91 30 i 

9140 Elcv cofftcts OISP * Elevation correction.* 

9130 FReKo “Tu*y*4Ptoiilf$tl,2J,Tdty(*) 

91d 0 FRUU *Tnit*EProi ilc$[l,2] ,Tnite(*) 

9170 FRCAD ”E lev* 4P rci } U$ 1 1 , 2 J , Topo (• ) 

91b0 J-Topc(l) 

9190 Dltt*Topo(2) 

920u L-0 

9210 INPol ‘Enter starting ano ending diitoncca for elcv cor rcl* «C1 « D2 

9220 11.D1/Dist4l 

9230 I2>02/Di&t.l 

9240 fOR 1-1142 TO 1242 

9230 L-L41 

9260 X(L)-Tc.ay(l) 

9270 y(L)-ropo(l) 

92dO Z(E) -iniied ) 

9290 x2(L)-{Tcay(l)4Tnite(l>)/2 
9300 Mat 1 

9310 KEOIM X(L) ,y(E) ,Z (E) ,X2 (L) 

9320 CALL La near ( V ( • ) , X {* ) , L , Al pna^oay ,Be ta_oay , De v_oay ) 

933u CALL Linear ( y (•), i (*)# 3,. ‘ lpna“ni te ,B e ta_ni te , Dev_^ni te) 

93a0 call Linear (y(*i ,X2 (*) ,L, a I pna mean, Beta mean, Dev mean) 

9350 printer is 0 
93t)0 FIXED 4 

9370 PRIM “Tuay- "eVALS (Alpna^day) a* • Elev 4 *i0AL3 (Beta_oay ) i* 4/- * 4 vALS ( Dtv^oay ) 

93bU PRINT "Tnite- ”4vALS(Aipna nite)4* • Elev 4 * 4vAL3 taela nite)4* */- *4VAL5(Dev nice) 

9390 PRINT "Tnean- *40AL3 (Alpna'/nc-an) 4* • Elev 4 *4VA1-5 (Beta]_mean ) 4* ♦/- * 4VAL5 (Dev_me£n) 

94 NJ PRINT 

9410 STANDARD 

94 20 PRIN'JLR IS 16 

943u PRINT L1N(3) 

9440 PRINT "Use CONT to rename TOay ano Tnit as Tot and Tnt * 

9450 PRINT "ano store elevation correction values in Toay and Tnit.” 

9460 PAUSE 

9470 RENAME *Toay*4P tof i leS 11 , 2) TO "Tdy * 4Prof i 1 eS ( 1 , 2 ) 

94d0 RENAME *Tni t" 4P rof i leS ( 1 , 2) TO *Tnt* 4Pt of i lej 1 1 , 2 ) 

9490 FCREATE *Tday ’ 4Pr of i le$ U , 2 ) , 52 
9500 FCREATE "Tni t" 4P rof 1 lej [1 , 2) , 52 

951i Elev_cori-l I Set elevation correction indicator 
•520 Alpna-Alpna mean 
9530 FOR 1-1 TO J 

9540 Taay{l42) -Toay(l42)-Alpna* (Topo(l42)-1400) 

•550 Tnite (142) -Tni te(l 42) -Alpna* (Topo{142)-l400) 

9560 NEXT 1 

• 570 FPRINT •Toay*4Proiile5U, 2J ,Tday(4) 

•Sao FPRINT *Tnit*4Prof ileS |1 , 2) ,Tnite (• ) 


PAGE IS 
CM'CV;' g:;M.lTY 


103 





IS 0 

tftwo nihT 

ftio SRINT *ro itstuit e»t» lilvfti* 
t420 PHINT “PUKCt Tu«y.. }R£NAI.L Tdy.. TO Td«y. . * 
.ts)d SRiin' *ruHue m^t.. iksname mt.. TOTnit.. * 
M40 I*R1NT 

nso rsihTCR is is 
fssd Disr "riKiSHCD.* 

SS70 GOTO Function keys 

fSiO I 

SS40 I 

9700 I 

9710 Rt^ort plot! I 

9720 DISP ^Report plotting.* 

9720 SERIAL 

9740 PLOTTER IS 7,S,*9S72A* 

97&0 COSUb .Mult^lilt read 

9760 MAT R«tl>Kct 1*(T00) i Convtrt back to t 
9770 D»>29 

9790 xii-*Distanct( km)* 

9790 K.AT SEARCH X(* ) ,HAX> Xmax 
9900 xminvO 

9610 .Oxa- (xmax-xmin) *. 126/6 
9920 DEG 

9930 /OR 1>1 TO 41 Pour antets 
9940 IF 1«1 THEN GOTO Plotfl 
9960 IF 1«2 THEN GOTO Plotk2 
9960 IF 1>3 THEN GOTO PlotaS 
. 9670 IF I«4 'rHEN GOTO Plott4 
9660 Plots, HAT V*Toay 
9690 y$»"Temp(K)* 

9900 Dy-6 

9910 Sp9»* - Day.* 

9920 HAT Duir^Tnitc 
9930 D5»"Temp(K)* 

9940 Doy-6 

9950 Sp2 6»* - Nignt." 

9960 GOIC 10230 

9970 Plots2: MAT Y«TopO 

9990 Y5»*Elevation(m) * 

999u Dy>100 

10000 Sp5-** 

10010 MAT Du.~«Detv 
10020 0$-*Gtradient (rn/km)* 

10030 Sp2»»"* 

10040 DoyalO 
10060 GOTO 10230 
lOUtO PlotsliKiAT Y-Dt 
10070 Y6«“Temp{K)* 

10090 Oy-6 

10090 Sp6»* - Temp oil I.* 

10100 KiAT Dun.>{4ean 
10110 D6-Y6 

10120 Sp2$»* - Temp mean.* 

10130 Day>Dy 
10140 GOTO 10230 
10160 Plots4:;4AT Y-Refl 
10160 Y$-*Re£l(l)* 

1017U Dy>6 
10160 Sp6-** 


/O'/ 



101 fU Ml bun«ln«fl 

1020U Duy«^0U 

1031U 0»*”lnvlil«(TlU)* 

1U320 Sp21-*" 

1021u rOH I'-l TO 2 
1U240 LUilT 0.4Ut,0,24S 
102^0 L.ORC & 

10240 CSUC 2.71 

1027b ir t>»i ThLU a1«(2*21.4-13.5)*u 0/3*l 
103»0 ir ><>2 7'liCN X)-Al«».l*21.4*)0i//«1 
1024U 0>X1«1.5*2».4*)OU/24& 

10300 ii>(2*2l.44)4)*)0O/3Xi> 

103)0 Y2>)r)40*21.4*10b/2X1 

103^0 nObt Xl4j*21.4*)bO/2*5,Y)-.5*21.>. *10u/24!> 

1033i> LAbCt X4 
1034U LOIfi »u 

1031O hOvt; Xl-.9*26.4*)yo/2»5,U44.2S»2S.4MOo/iil 
lU3bu U THLtj Y4«D4 
10 370 UL4L VI 
lU3ou LOIF 0 

lu3yu XI ,X2,Y) ,Y2 

104UO rKA<’>&. 

104 lo LCiRG S 

1O420 IF »»»2 IritN I'AT Y*Dur. 

Io430 IF F*2 ThCI« Cy*Doy 
10440 If P-2 'intN lif.>*i>{24 
10«3U I.A'i' StAHCri Y(*) |NAX>Y(»4X 
lU4i)0 i.Ai' btMACri Y (* ) .KIN; Yn.i n 
10470 Ym n-2M ( Yrtiin/Dy ) *Dy 
104o0 Y«r.««>(2N'i (Y>['««/Oy) 4))*Dy 
1 u4 yo D/S - ( Yri«*-Sn\i n) * . 12 1/X . 3 
10 loo SCALL An.in,Am«»,yirin,y.T)£.x 
luDO AALS £>«,Dy, Ajnin,Yn.in.4,2 
10S20 OolPU 705 j"2S2>" 

10130 hOOL XininiYn-.ax 
10140 CRA«< Xn.aA,Ym«« 
lOSlO DRAW XraXtYmin 

lOloo FOR lc-1 TO 0 I Draw curv# 

10170 IF Ic-1 THEN FOvE X(IC),Y{lC^ 
lOlaO DRAn X(lC) . Y(lc) 
lOSlU NEXT IC 
lOaoO OOlPur 705;"VNj* 

lOtIO CSIZE 2 I LaDtl axes 

10120 LORG S 

10b30 FOR ya-ymin TO Ymax STEP Dy 

10o40 F.OVE X<‘nrn-Dxa, Ys 

10610 LABEL WAL$(Ys) 

lu»60 NEXT Ya 

10670 LORG 6 

lOioO FOR Xa-Xmin TO Xmax STEP Dx 
lObtfu NOVL xa.Ynirt-DyK 
107u0 LABEL VALl(Xk) 

10710 NEXT X» 

10720 LORG 5 

10730 KOVE (xnln4Xmax)*.l,Yn^x>Dy»*2 

10740 LABEL •Pioxxia • XProl i leJll , 2] 6* 6Sp5i* •6Area$6“ Area* 

10 710 PEN 0 

10760 NEXT P 

10770 OUTPUT 70lj*AHr* 

10760 OUTPUT /03j*AH;‘ 


/O^ 


107f;V hCXT I 

)U»UO OUTi*UT 706(*Ant* 

)0§Jb GOTO f'unctton 
ifliSi/ I 

.]0»JO I 

I0«40 I 

lOk&U ttUM Linc«MA(* ) , V (•) ,L.Alpn«,ttl«,Dtv) 
lUkvU I ril V • Alpn**A 4 ItlA 4/- Dtv 
J0»7u OPTIO.. MASC 1 

]0*«U D1.4 A) (5J0 J) ,yj (*»^i<.l),Atrns() ,S(/U) .0( i .) ) . r^( 1 • 1 ) .r (1 , ) ) 

)b»»u AtOlM Aiail) .yj (L,J),Atrnf(l>U 
J0Vw(/ rOH l.J 10 L 
ioyjo Ai u.i).A(i) 

J0(2i) yj (1 JJ>M1 ) 

JUfcJi. l.lAl 1 

J094l« An«dn.SU.4( K!) 

)oysi) ymt«(i.sJK(y))/L 
lUvCb hAT Al.Al-iA.'ntAn) 

10»7t KA1 yl.y]>n.Tt«n) 

JuybO f.AT Aun».r>itf (Ai ) 
luyyw Mi D.Atrni*xi 
llOtUi Mr A>ltiv(0) 
iiwio Mr o.Atins*yi 
M1‘ r.A'O 
llOJij AlpriA-r (1 « J) 

))04U 0c iA.yr«*n-AlpDA* Air.-4n 
J)06u iiAl a 1 .Al4 (a.?caii) 
f-AT yj.y J4 (Vme*n) 
iJ07C I:mT aJ-aJ* (Al^ni) 

Mr Al.Al4(0ttt) 
lioyo MAT Al-Al-yl 
llluu Ml Al>Al .A] 
nuo DCV>S.N(SuiMAl)/(t-))) 

)112U SObLNO 

lllji) SUi> rAAt liiit (X(*) fOkiL.Aupper ,Alo%ier,A.TiR,ArAA) 
iri.O I CutpUc'iA AA>ll<,AnA« 

11)^0 OPTION oASL 1 

IHou K£OIil A(LJ 

J1J7U A..au>:>Utl( A)/L. OlV OA*OA 

lllev) i-.Ar bCAACH A(* ) «LOC ( >-Auppvr> ,’NuppCI 

1J1»U A.»AA.An<iu40x 

112UU MAX »*Atcn-.MAT SLARCH A (* } «LOC (>XmAX) {N 
ir A-Nupi:«r<L*.0< TnCU run 
1122u A.TiAX .AtnAx 40x 
l)23vi GOTO MAX BCArcn 

JJ240 MinsMAT iLARCrt A ( • ) ,LOC (,<« Alowcr ) i W lower 
1125U ARU n- Airin-Ox 

llTbO Hin^tcArcniAiAT SCARCri X(*) ,LOC(<A.'nin) >N 

1127U IF A-Nlowcr <L* .U4 THEN Finiin 

112a0 AAixn.XTiin-Dx 

l)2y0 GOTO Min^BCAicii 

llJOO Pi niAnssUBCNO 

11310 END 
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